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(57) Abstract: NMR methods for extracting information about a fluid in rock and logging apparatuses for implementing such meth- 
ods in a borehole environment are provided. The methods involve generating at least two different magnetic field pulse sequences. 
The magnetic field pulse sequences include a first portion and a second portion. A magnetic field pulse sequence is generated, and 
magnetic resonance signals are detected using the second portion of the sequence. The first portion of the sequence is modified, and 
again the sequence generated and magnetic resonance signals detected using the second portion. The magnetic resonance signals 
are analyzed, and information about, for example, diffusion coefficient, viscosity, composition, saturation in a rock, pore size, pore 
geometry and the like, extracted from the analyzed signals. 
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NUCLEAR MAGNETIC RESONANCE METHOD AND LOGGING APPARATUS 

This patent application claims priority from U.S. Provisional Application No. 
60/170,121 filed on December 10, 1999 and from U.S. Application Serial No. 
5 09/528,881 filed on March 20, 2000, both of which are herein incorporated by 
reference in their entireties. 

FIELD OF THE INVENTION 
This invention relates to investigations of rock samples, and more particularly 
10 relates to nuclear magnetic resonance (NMR) methods for determining characteristics 
of subsurface rock, including fluid composition. 

BACKGROUND 

NMR has been a common laboratory technique for over forty years and has 
15 become an important tool in formation evaluation. General background of NMR well 
logging can be found, for example, in U.S. Patent 5,023,551 to Kleinberg et aL, which 
is assigned to the same assignee as the present invention and herein incorporated by 
reference in its entirety. 

NMR relies upon the fact that the nuclei of many chemical elements have 
20 angular momentum ("spin") and a magnetic moment. In an externally applied static 
magnetic field, the spins of nuclei align themselves along the direction of the static 
field. This equilibrium situation can be disturbed by a pulse of an oscillating 
magnetic field (e.g., an RF pulse) that tips the spins away from the static field 
direction. The angle through which the spins are tipped is given by 9=YB 3 tp/2, where 
25 y is the gyromagnetic ratio, B] is the linearly polarized oscillating field strength, and t p 
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is the duration of the pulse. Tipping pulses of ninety and one hundred eighty degrees 
are most common. 

After tipping, two things occur simultaneously. First, the spins precess around 
the direction of the static field at the Larmor frequency, given by coo=YB 0 , where B 0 is 
5 the strength of the static field and y is the gyromagnetic ratio. For hydrogen nuclei, 
Y/2tc=4258 Hz/Gauss, so, for example, in a static field of 235 Gauss, the hydrogen 
spins would precess at a frequency of ] MHz. Second, the spins return to the 
equilibrium direction according to a decay time, Ti, which is known as the spin-lattice 
relaxation time. 

10 Also associated with the spin of molecular nuclei is a second relaxation time, 

T 2 , called the spin-spin relaxation time. At the end of a ninety-degree tipping pulse, 
all the spins are pointed in a common direction perpendicular, or transverse, to the 
static field, and they all precess at the Larmor frequency. However, because of small 
fluctuations in the static field induced by other spins or paramagnetic impurities, the 

15 spins precess at slightly different frequencies, and the transverse magnetization 
dephases with a time constant T 2 - 

A standard technique for measuring T 2 , both in the laboratory and in well 
logging, uses an RF pulse sequence known as the CPMG (Carr-Purcell-Meiboom- 
Gill) sequence. As is well known, after a wait time that precedes each pulse 

20 sequence, a ninety degree pulse causes the spins to start precessing. Then, a one 

hundred eighty degree pulse is applied that keeps the spins in the measurement plane, 
but causes the spins, which are dephasing in the transverse plane, to reverse direction 
and to refocus. By repeatedly reversing the spins using a series of one hundred eighty 
degree pulses, a series of ''spin echoes" appear. The train of echoes is measured and 

25 processed to determine the irreversible dephasing, T 2 . 
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In rock formations, such as in a borehole environment, T 2 for hydrogen- 
containing fluids (e.g., water, oil, gas) can have significant contributions due to 
surface relaxation, bulk relaxation, and diffusion effects, i.e., 

1111 
— = + + (1) 

J 2 1 2, surf ace 1 2,bulk 1 2 .diffusion 

5 Each of these contributions provides some information about the rock formation 
and/or about the fluid in the rock formation. For example, in a wetting phase, the 
surface relaxation contribution, T2,surface, dominates the distribution of observed 
distribution of decay times,/^). Spins relax predominantly due to collisions with 
the grain surface, with the collision rate being inversely proportional to the pore size. 

10 This means that the observed relaxation time is roughly proportional to the pore size, 
i.e., lAr2 iS urface = teS/Vp, where S is the surface area of the pore, V p is the pore volume, 
and p 2 is the surface relaxivity of the rock, a phenomenological parameter that 
indicates how relaxing the surface is. Thus, for a wetting phase, the observed f(J 2 ) 
provides information about pore size distribution. In a nonwetting phase, surface 

15 relaxation becomes negligible and bulk relaxation, which is related to viscosity, 
dominates the observed j^T 2 ). Thus, for a nonwetting phase, the observed flj 2 ) 
provides information about viscosity. 

In a uniform static magnetic field, each spin will experience the same 
magnetic field strength regardless of its position within the static field, and diffusion 

20 will not contribute to the observed flj 2 ). In a magnetic field gradient, however, each 
spin will experience different magnetic field strengths as it diffuses through the static 
field. The Larmor frequencies of the diffusing spins become time dependent, and the 
series of one hundred eighty degree pulses cannot refocus the spins completely, 
leading to an additional decay signal. This additional decay signal is proportional to 
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the diffusion coefficient, D, of the fluid and to the square of the gradient strength, g, 
and the square of the echo spacing, t£> i.e., 



As the diffusion coefficient provides an indication of fluid type, measurement of the 
5 diffusion effects on/Tr 2 ) can be used as the basis for determining the types of fluids in 
a rock formation. 

Certain NMR measurements of diffusion involve changing the echo spacing, 
t E , in a standard CPMG sequence, and thus the amount of diffusion the spins undergo 
between echoes, and then comparing the measured relaxations. Figures 1 A and IB 

10 generally illustrate this approach. Figure 1A shows two CPMG sequences with 
different echo spacings, ti and t 2 , where t 2 is longer than ti. As the echo spacing 
increases, the spins diffuse further between echoes, and the measured relaxation times 
will decrease depending on the diffusion coefficient of the fluid, as given in Equation 
2 above. Figure IB shows the relaxation distributions, flT 2 ), for an oil and water 

15 determined from the two sets of echoes acquired from the two CPMG sequences 
illustrated in Figure 1A. As seen in Figure IB, the relaxation distribution with the 
longer echo spacing, t 2 , is shifted to lower relaxation times, T 2 , relative to the 
relaxation distribution with the shorter echo spacing, ti. The size of the shift is 
proportional to the size of the diffusion coefficient, as indicated by arrows 1 and 2. 

20 The shift of ftT 2 ) for a fluid with a small diffusion coefficient 1, such as heavy oil, is 
smaller than the shift for a fluid with a larger diffusion coefficient 2, such as water or 
natural gas. 

While such NMR diffusion measurements can be useful, they suffer from a 
number of drawbacks. For example, for a given acquisition time, the two CPMG 



2, diffusion 




(2). 
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sequences will not have the same number of echoes. The CPMG sequence with 
longer echo spacing will have a fewer number of echoes available, so will suffer from 
lower signal to noise and lower data quality in general. In addition, relaxation 
distributions for different fluids often overlap, at least partially, making it difficult to 
identify shifts of individual relaxation times. In cases where the diffusion coefficients 
for different fluids are small, the shifts may be difficult to distinguish. Finally, these 
methods cannot separate out the contributions due to diffusion effects from the 
surface and bulk relaxation contributions in the observed relaxation distributions. 
Surface relaxation and diffusion have similar effects on the observed relaxation 
distributions, so these methods may provide inaccurate information about the fluid 
and about the rock or earth formation under investigation. 

SUMMARY OF INVENTION 

The invention provides in one aspect a method for extracting information 
about a fluid that may be contained in rock or within a portion of earth formation 
surrounding a borehole (as used hereinafter, the term "rock" includes earth, earth 
formation, and a portion of earth formation). For example, certain embodiments of 
the invention involve generating a sequence of magnetic field pulses in the fluid. The 
magnetic field pulse sequence includes an initial magnetic field pulse, a first portion 
that follows the initial magnetic field pulse, and a second portion that follows the first 
portion. Magnetic resonance signals are detected using the second portion. The first 
portion of the sequence is then modified, and the magnetic field pulse sequence 
generated and magnetic resonance signals detected using the second portion again. 
The magnetic resonance signals are analyzed from a time relative to the initial 
magnetic field pulse, and information about the fluid extracted. 



WO 01/42817 ^pT/USOO/32589 

6 

A second aspect of the invention provides a logging apparatus for 
implementing embodiments of the inventive measurements in a borehole 
environment. One embodiment of such a logging apparatus includes a logging tool 
that is moveable through a borehole and a processor coupled with the logging tool. 
5 The processor is programmed with instructions which, when executed by the 

processor, cause the logging tool to generate a sequence of magnetic field pulses in a 
region of investigation of earth formation surrounding the borehole. The magnetic 
field pulse sequence includes an initial magnetic field pulse, a first portion, and a 
second portion. The processor causes the logging tool to detect magnetic resonance 
10 signals produced from the region of investigation using the second portion of the 

sequence, and then modify the first portion of the sequence and repeat generating the 
sequence and detecting magnetic resonance signals. The programmed instructions 
also cause the processor to analyze magnetic resonance signals from a time relative to 
the initial magnetic field pulse and extract information about the region of 
15 investigation. 

Further details, features and embodiments of the invention will become more 
readily apparent from the detailed description that follows. 

BRIEF DESCRIPTION OF FIGURES 
The invention will be described in more detail below in conjunction with the 

20 following Figures, in which: 

Figures 1A and IB, taken together, illustrate a NMR measurement and T 2 

distributions obtained therefrom according to the prior art; 

Figures 2 A and 2B, taken together, illustrate one embodiment of a NMR 

measurement and T 2 distributions obtained therefrom according to the invention; 
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Figure 3 shows a modified CPMG sequence that can be used in a NMR 
measurement according to embodiments of the invention; 

Figure 4 is a graph of attenuation versus initial echo spacing for three different 
fluids according to certain embodiments of the invention; 
5 Figure 5 is a graph showing relative signal attenuations for three samples as 

initial echo spacing is varied according certain embodiments of the invention; 
% Figure 6 illustrates a plurality of magnetic field pulse sequences that can be 

used in one embodiment of a NMR measurement according to the invention; 

Figure 7 shows two-dimensional plot of diffusion coefficient versus relaxation 
1 0 time according to certain embodiments of the invention; 

Figure 8 illustrates a plurality of magnetic field pulse sequences that can be 
used in another embodiment of a NMR measurement according to the invention; 

Figure 9 shows two graphs illustrating results from two different rock samples 
of a time dependent diffusion coefficient measurement according to certain 
15 embodiments of the invention; 

Figure 10 is a schematic diagram, partially in block form, of one embodiment 
of a well logging apparatus that can be used in implementing methods according to 
the invention in a borehole environment; 

Figure 1 1 is a block diagram of downhole circuitry that can be used with a 
20 well logging apparatus, such as shown in Figure 10; and 

Figure 12 is a flow diagram of a routine that can be used in programming a 
. processor or processors in implementing certain embodiments of the invention. 
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DETAILED DESCRIPTION 

A NMR measurement according to the invention involves generating at least 
two different magnetic field pulse sequences hi a fluid and detecting magnetic 
5 resonance signals produced by the different magnetic field pulse sequences. 

Information, such as saturation, diffusion coefficient, viscosity, composition, etc., 
about a fluid, e.g., in a rock sample or in a rock formation surrounding a borehole, is 
extracted by analyzing the different magnetic resonance signals. 

Generally speaking, the magnetic field pulse sequences used in the invention 

10 can be characterized as having two portions, a first portion followed by a second 

portion. In a NMR measurement according to the invention, the sensitivity of the first 
portion to diffusion effects in the presence of a magnetic field gradient is modified 
while the second portion remains substantially the same. Magnetic resonance signals 
are detected using the second portion and analyzed. Signals detected using the second 

15 portion, which is not changed, will exhibit substantially the same relaxation 

distribution from one magnetic field pulse sequence to the next, except that the 
amplitude of the signals will have been altered as a result of the first portion having 
been modified. By analyzing how the magnetic resonance signals change as the first 
portion is modified, information about the fluid may be extracted. 

20 Stated another way, in the presence of a magnetic field gradient, the first 

portion of the magnetic field pulse sequences used in the invention are sensitive to 
both relaxation and diffusion effects whereas the second portion has substantially the 
same sensitivity to relaxation effects but reduced sensitivity to diffusion effects. 
Diffusion effects during the first portion introduce an extra decay into the signal 

25 detected using the second portion. The signal detected using the second portion is 
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thus attenuated, or "diffusion edited" in proportion to die diffusion coefficient of the 
fluid (see Equation 2, above). 

These diffusion effects may be detected in the presence of a static magnetic 
field gradient, or with the use of pulsed field gradients as described, for example, in 
5 U.S. Patent No. 5,796,252 to KJeinberg et aL, which is assigned to the same assignee 
as the present invention and incorporated herein by reference in its entirety. Pulsed 
^ field gradients introduced into the first portion of the magnetic field pulse sequences 
^ of the invention also may be used in conjunction with a static magnetic field to 
enhance these diffusion effects. 
10 Figures 2 A and 2B illustrate one embodiment of a NMR measurement 

according to the invention. After a static magnetic field is generated in a fluid in a 
rock, a first magnetic field pulse sequence 3 and then a second magnetic field pulse 
sequence 4 are generated in the fluid. The first sequence 3 in this embodiment 
(shown at the top of Figure 2A) is a standard CPMG sequence in which an initial 90- 
15 degree pulse is followed, after a time t^n, by a series of 180-degree pulses separated 
in time by about 2t sh0 ji- A magnetic resonance spin echo appears a time t sh0 rt after 
each 1 80-degree pulse, producing a series of magnetic spin echoes with a time 
spacing, t E ,sb 0 n, approximately equal to 2^^. This first sequence may be represented 
as: 

20 90 - [t^ - 1 80 - t shon - echoj n (3) 

where the time separating the echoes, t E .d,an, is equal to about 2^; echo; is the i* 
magnetic resonance spin echo; and n is the number of spin echoes in the sequence. 

A standard CPMG sequence may be characterized as a magnetic field pulse 
sequence according to the invention in which the first portion is substantially identical 

25 to the second portion. A modified CPMG sequence according to the invention may 
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be thought of as a CPMG sequence in which the first portion has been modified. The 
second portion of the modified CPMG sequence is not changed and so generates a 
CPMG-like series of magnetic resonance spin echoes with a time spacing 
approximately equal to tE.shon- 

The second sequence 4 (shown at the bottom of Figure 2A) is an embodiment 
of a modified CPMG sequence 4 in which the first few echo spacings of the standard 
CPMG sequence are elongated. An initial 90-degree pulse is followed by a first 
portion containing a first series of 180-degree pulses that begins a time, ti OJ i g .i, after the 
initial 90-degree pulse and are separated by about 2tj OJ ig.i, where t) ong> i is greater than 
tshon. A magnetic resonance spin echo appears at a time t] on g,i after each 180-degree 
pulse, producing a first series of magnetic resonance spin echoes with a time spacing 
t E ,i approximately equal to 2t] on o,j. The first portion is followed a second portion 
containing a second series of 180-degree pulses separated in time by about 2tshort- The 
second series of 180-degree pulses begins at a time t^on after the last spin echo of the 
first portion and refocuses this last spin echo to produce a second series of magnetic 
resonance spin echoes having a time spacing, t E ,shojt, which is equal to about 2tshort. 

The embodiment of a modified CPMG sequence 4 shown at the bottom of 
Figure 2A may be represented generally as: 



where, for the j sequence, t] 0i i g j is greater than t^on, echo^j is the k magnetic 
resonance spin echo of the first portion; (m,j) is the number of spin echoes of the first 
portion; echojj is the i ,th magnetic resonance spin echo of the second portion; and 
(n\j) is the number of spin echoes of the second portion. As shown in Figure 2A, the 
first portion of the modified CPMG sequence 4 contains two spin echoes, i.e., 



90 - [tjong j - 1 80 - tjongj - echOj: j] m j - [tshort -180- tshon - 



echorjL-j (4) 
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(m,j)=2; it will be appreciated, however, that the first portion may have other numbers 
of echoes. 

Another embodiment of a modified CPMG sequence that may be used in 
accordance with the invention has a first portion that contains a stimulated echo 
5 sequence. For example, as shown in Figure 3, the first portion includes a series of 
two 90-degree pulses at times 5j and Ai after an initial 90-degree pulse, producing a 
stimulated echo at a time (A^S^ after the initial 90-degree pulse. The second 
portion, which follows the first portion, contains a series of 180-degree pulses which 
begins a time tshon after the last stimulated echo, refocusing the stimulated echo to 
10 produce a series of magnetic resonance spin echoes having a time spacing, fe^on, 

which is equal to about 2tsh on . The overall sequence may be represented generally as: 
90 - [5j - 90 - (A r 5j) - 90 - Sj -echo k jj^j - [t^ - 1 80 - Wt - echorjVjCS) 
where, for the j* sequence, echo^j is the k* stimulated echo of the first portion; (m,j) 
is the number of stimulated echoes of the first portion having an echo spacing equal to 
15 about (Aj+5j); and 8j, Aj, tshon, and indices and (n\j) are as defined above. 

Other embodiments of a modified CPMG sequence that may be used in 
accordance with the invention are described in U.S. Application Serial No. 
09/528,881, filed on March 20, 2000, which is incorporated herein by reference. 

The graph in Figure 2B represents the relaxation distributions extracted from 
20 the magnetic resonance signals detected using the second portions of the first 

sequence 3 and the second sequence 4 in the presence of a magnetic field gradient. 
The relaxation distributions,^^), for an oil and water, shown in Figure 2B are 
identical except that the amplitude of the second sequence signal 4 is attenuated 
relative to the amplitude of the first sequence signal 3. The amount of the signal 
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attenuation is proportional to the size of the diffusion coefficient, as indicated by 
arrows 5 and 6. Thus, the signal attenuation for a fluid with a small diffusion 
coefficient 5, such as oil, is smaller than the signal attenuation for a fluid with a larger 
diffusion coefficient 6, such as water or natural gas. Measuring the relative change in 
5 amplitude between corresponding spin echoes from the first sequence and from the 
second sequence can yield quantitative information about the fluid in the sample 
under investigation. 

Typically, the echo spacing in the second portion, tE tS hoit, is chosen to be as 
short as possible to increase the number, n, of spin echoes that can be generated and 

10 detected within a given acquisition time. This, in general, increases the signal-to- 
noise of the measured signal and, in the presence of a magnetic field gradient, reduces 
the sensitivity of the second portion to diffusion effects. Times for tsh ort on the order 
of about 0. 1 milliseconds (100 jxs), leading to echo spacings, t^shoit, on the order of 
about 0.2 ms (200 jxs), are currently used in well logging measurements, although the 

15 measurements of the invention may be made using other times for tshon and t E ,shoit- 
To properly correlate the spin echoes from the first sequence with the spin 
echoes from the second sequence, the data processing for both sequences starts at the 
same time, ta, relative to the initial 90-degree pulse. In other words, spin echoes that 
precede t<j are not used in processing, and only spin echoes starting from t<j are 

20 analyzed. By beginning the data processing at the same time for each sequence, the 
magnetic resonance signals from each sequence will reflect substantially similar 
relaxation effects. 

As shown in Figure 2A, U corresponds to the last echo of the first portion of 
the modified CPMG sequence. For example, with tshojt^O.l ms and t] OJlg j=4 ms, data 
25 processing would start at t s =16ms, and if the two spin echoes of the first portion of 
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the second sequence are disregarded, then the first 80 spin echoes of the first 
sequence, Le., echo] from i=l to (m,j)x(ti ong j/t sho n), would be disregarded. Data 
processing according to the certain embodiments of the invention, however, may 
include the last echo of the first portion, i.e>, the echo from the first portion that is 
5 refocused by the magnetic pulses of the second portion. This means that the spin 
echoes that occur at t=t<j would be included in the data processing. 

Both sequences also may be truncated at the end during processing so a 
? substantially equal number of spin echoes are analyzed and/or echoes with poor 
signal-to-noise are disregarded. Typically, both sequences still will contain a large 
10 number of echoes, particularly as compared to prior art techniques in which the echo 
spacing of a standard CPMG is increased. Thus, the methods of the invention will 
generally allow more precise measurements and provide better signal-to-noise as 
compared to prior art techniques. 

Though the spin echoes preceding U typically are not used in the above 
15 analysis, they may be collected and analyzed to extract other information, such as 
total porosity, permeability, etc. of a fluid-containing rock sample, as described, for 
example, in U.S. Patent Nos. 5,023,551, 5,363,041, and 5,389,877. Such analyses 
may be done in conjunction with the methods of this invention. 

A diffusion edited signal generated according to (4), for times > t<j, may be 
20 represented as: 

M(t EJ ,t) = JJdDrfr 2 /(Ar 2 )exp(-^)exp(-^ r ^ 2 Z>4..) (6), 

where u coiresponds to the time (relative to the initial pulse) of the last echo of the 
first portion, y(D,T 2 ) is the two-dimensional diffusion-T 2 probability density function, 
(m,j) is the number of echoes of the first portion, y is the gyromagnetic ratio, g is the 
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gradient field strength, D is the diffusion coefficient, and t£j is the echo spacing in 
the first portion of the modified sequence. For modified sequences having more than 
one echo in the first portion, Le., (m,j)>l, multiexponential diffusion decays may be 
observed over a broad detection bandwidth. Such multiexponential decay can be 
5 modeled more accurately by replacing the single exponential diffusion attenuation 
term in Equation 6 with multiple terms. For example, for (mj)=2, two exponential 
terms, one from the direct echo and the other from the stimulated echo, would replace 
the single exponential term in Equation 6, yielding: 

M(/ £J ,o = JJ^r 2 ^ (7), 

10 where A and B are parameters that depend only on the detection bandwidth. 

Similarly, a diffusion edited signal generated according to (5), for times > t<j, 
may be represented as: 

M u i5j) = ±fjdDdT t f(D0^T 2 )^ (8), 

2 ^2 

where t<j is the time (relative to the initial pulse) of the stimulated echo.,y(D(td),T2) is 
15 the two-dimensional diffusion-T 2 probability density function at time t<j, and D(ta) 
represents the diffusion coefficient at time t<i. 

The amplitude of the detected echoes can be measured using any of various 
signal processing techniques known in the art and then, according to certain 
embodiments of the invention, fit to one of Equation 6, 7, 8, or other diffusion edited 
20 signal equations, depending on the magnetic pulse sequence used, in order to extract 
diffusion coefficient and other information about the fluid. 

As previously suggested, measuring the relative change in amplitude between 
corresponding spin echoes from the first and second magnetic field pulse sequences 
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can provide quantitative information about the fluids in the sample under 
investigation. In particular, a comparison of corresponding spin echoes can separate 
different contributions to the T 2 relaxation distribution, and so yield more accurate 
. information about the fluid. For example, as discussed* with regards to Figure 2A, in 
5 the presence of a magnetic field gradient, the first portion of the modified CPMG 
sequences of the invention are more sensitive to diffusion effects than the second 
portion, leading to the signal from the second portion to be attenuated, or diffusion- 
edited, compared to the standard CPMG signal. However, the signal from the second 
portion exhibits the same surface and bulk relaxation effects as, and is otherwise 
10 identical to, the signal from the standard CPMG sequence. Thus, taking a ratio of the 
amplitudes of the corresponding T2 distributions of the first and second sequences 
separates diffusion effects from both surface and buUc relaxation effects. 

For multiphase fluids, linear mixing laws govern the observed signal 
attenuation. Thus, components of the fluid having different diffusion coefficients 
15 may be differentiated by appropriately selecting and varying a time spacing in the first 
portion of the modified sequences (e.g., t E ,i; 8j or Aj). Figure 4 contains a graph that 
may be useful in selecting a first portion time spacing. The graph shows approximate 
attenuation versus initial echo spacing, t E ,i, for oil (having a diffusion coefficient, D^, 
of 2.3 x 1(T 6 cm 2 /s), water (D watM =2.3 x 10" 5 cm 2 /s) and gas (D gas =1.2 x 10' 3 cm 2 /s) 
20 measured in a magnetic field gradient of about 25 Gauss/cm with a modified sequence 
generated according to (4), having two echoes in the first portion. Over a range of 
echo spacings greater than approximately 0.8 ms, the relative attenuation of the gas 
signal is much larger than for either water or oil, indicating that this range may be 
useful in differentiating gas from water or oil. In general, an initial echo spacing in a 
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range of between about 0.5 to about 20 ms may be useful for hydrocarbon typing 
applications of the invention. 

An initial echo spacing to differentiate between two diffusing components 
may be approximated as: 



D / nl/3 



y 2 g 2 (D + -D_) 



(9), 



where D+ is the diffusion coefficient of the higher diffusing component, D. is the 
diffusion coefficient of the lower diffusing component, y is the gyromagnetic ratio and 
g is the magnetic field gradient strength. Equation 9 may provide a useful starting 
point in selecting an initial time spacing, if the different diffusing components (or a 

10 range of diffusion coefficients) are known or can be approximated a priori. 

In some embodiments of the invention, the relative attenuation of the echo 
signal may be quantified in terms of an attenuation factor, ai. The fluid composition 
then may be determined by taking a linear combination of the attenuation factors for 
the individual components. For example, the attenuation factor may be calculated as 

15 a ratio of the sums of all the measured echoes of a modified CPMG sequence to the 

measured echoes of a standard CPMG sequence, or as a ratio of initial amplitudes of a 
modified CPMG sequence to a standard CPMG sequence. In a gradient field of about 
25 Gauss/cm and using a modified sequence generated according to (4) with two 
echoes in the first portion (t E ,i = 8 ms, t E ,shojt= 0.2 ms), an attenuation factor for bulk 

20 water, a w , was calculated to be about 0.32 and for 6 cp oil, was calculated to be 

about 0.85. For a fluid having water and oil components, the water saturation, S w , i.e., 
the proportion of fluid that is water, can be determined from the measured attenuation 
factor of the fluid, a meas , using the following relationship: 
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S w = a ° u ameas (10). 
a , — a 

oil w 

According to other embodiments, the fluid composition can be determined by 
fitting the measured attenuation curve versus a first portion time spacing as a 
superposition of component curves. Figure 5 shows the relative initial attenuations of 
5 diffusion edited signals as an initial echo spacing is changed for a NiCl-doped water 
sample (labeled "water"), a S6 oil sample (labeled "oil"), and a partially oil saturated 
sandstone sample. The diffusion edited signals were generated in a uniform gradient 
field of about 13 Gauss/cm with a modified CPMG sequence according to (4) and 
having two echoes in the first portion. For the water and oil samples, the plotted 
10 points, representing the amplitude of a modified CPMG signal relative to a standard 
CPMG signal, were fit to a curve using Equation 7 above, with A=0.59 and B=0.20. 
From the fit, the diffusion coefficients were extracted and found to be D v .=- XlO 5 
cm 2 /s for the water sample and D o ij=1.35xl0" 6 cm 2 /s for the oil sample. Tot he 
partially oil saturated sandstone sample, the relative attenuation was modeled as a 
15 superposition of two curves using Equation 7 and three parameters: S w , the n ater 

saturation; D w>e ff, the effective diffusion coefficient for water; and D 0 ii, e ff. the i ffective 
diffusion coefficient for oil. From the fit, it was determined that: S w =0.21, 
D w ,eff=1.9xl0" 5 cm 2 /s, and Dou.eff=13xl0" 6 cm 2 /s. 

The water saturation, S w , determined in this manner agreed with independent 
20 gravimetrical measurements, and the effective oil diffusion coefficient D , c matched 
the bulk diffusion coefficient D 0 a of the pure fluid, as expected. The effective water 
diffusion coefficient D w , c ff, however, was smaller than the bulk diffusion coe. ricient, 
D w , of water. This difference in water diffusion coefficients is attributed to restricted 
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diffusion of tlie water in the rock, as water in sandstone is known to mainly occupy 
the small pores of the rock. 

The effects of pore geometry, in terms of restricted diffusion and internal 
gradients, may be addressed using the methods described in U.S. Application Serial 
5 No. 09/528,881 . The effects of restricted diffusion and internal gradients may also be 
addressed in terms of the measured attenuation factor. A hypothetical water 
saturation for water that does not contribute to the observed signal attenuation, 

calculated as, e.g., S w res1r = a ™<* ~° w , 

and a ratio of the measured attenuation factor to that of bulk water, i.e., ameas , may be 

10 used to correct for the effects of restricted diffusion and internal gradients, 

respectively, in calculating, for example, water/oil saturation in a rock sample. 

The diffusion editing embodiments of the invention thus can extract saturation 
and diffusion coefficient information about a fluid directly without a priori knowledge 
of, and without having to make any assumptions about, relaxation in the fluid. As the 

15 measurements of these embodiments take place in the presence of a magnetic field 
gradient, the spin echoes are generated from a thin slice across the sample. By 
moving the sample and the gradient relative to one another along the gradient 
direction, an attribute profile of the fluid, for example, a saturation profile of a fluid in 
a rock sample, can be obtained. 
20 In other embodiments of the invention, a plurality of magnetic field pulse 

sequences, each with a different first portion (e.g., different time spacing, ti OJlgJ ; or 5j 
and/or Aj), are used in making the NMR measurements. The diffusion edited 
magnetic resonance signal depends on a time spacing in the first portion (e.g., t£j; Sj, 
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Aj), so fluids with different diffusion coefficients will diffuse different amounts during 
different time spacings (see Equations 6, 7, or 8, above). Inverting data using only 
two different time spacings (as described in the embodiments above) yields an 
average diffusion coefficient for every relaxation time and so may not resolve 
5 different diffusing components that have overlapping relaxation times. Measuring 

magnetic resonance signals using more than two different time spacings allows 
* diffusion coefficient and relaxation time to be extracted separately and so can help 
resolve different fluid components in terms of diffusion and relaxation, as well as 
other parameters that may be derived from D, T 2 , or a combination of D and T 2 , such 
10 as Ti, viscosity, saturation, porosity, etc.. 

Figure 6 illustrates a plurality of magnetic field pulse sequences that may be 
used in one embodiment of a NMR measurement according to the invention. A first 
magnetic field pulse sequence 10 is a standard CPMG sequence having a first time 
spacing, t^*, and a first echo spacing, t EfS hort. A second magnetic field pulse sequence 
15 11 is generated by modifying, in this case elongating, a time spacing between the 
initial 90-degree pulse and the first 180-degree pulse, t l0 j, g ,i, leading to an elongated 
echo spacing, t E ,i, in the first portion. A third magnetic field pulse sequence 12 is 
generated by modifying a time spacing of the first portion again, tj OJlg , 2 . A fourth 
magnetic field pulse sequence 13 is again generated by modifying a time spacing of 
20 the first portion, ti png 3> and also the number of spin echoes in the first portion. 

This embodiment of a NMR measurement involves generating each of the 
magnetic field pulse sequences shown in Figure 6 in a fluid in a rock. Magnetic 
resonance signals from each sequence are detected using the second portion of each 
sequence, which remains substantially the same from sequence to sequence.. As 
25 discussed above, to properly correlate the spin echoes detected from each sequence, 
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the data processing for each sequence starts at a time U relative to the initial 90-degree 
pulse. Figure 6 shows U corresponding to the time of the last echo of the second 
sequence 11, which corresponds in this illustration with the time of the last echo of 
the fourth sequence 13. Some spin echoes from the second portion of the third 
5 sequence 12, which arise prior to time t<j, are discarded in the data processing 

according to this embodiment. The time td need not necessarily correspond with the 
time of the last echo of the longest first portion, as shown in Figure 6; however, td 
typically is at least as long as the longest first portion. 

Although Figure 6 shows four magnetic field pulse sequences, other numbers 

10 of magnetic field pulse sequences may be used in NMR measurements according to 
the invention, with the use of more sequences resulting in higher resolution in the 
diffusion (or relaxation, viscosity, saturation, etc.) distribution. Additionally, the 
plurality of magnetic field pulse sequences need not include one type of modified 
CPMG sequence, as shown in Figure 6. Measurements according to methods of the 

15 invention may be made using a plurality of magnetic field pulse sequences that 

includes a combination of modified CPMG sequences generated according to (4), (5) 
or any of the other sequences described in U.S. Application Serial No. 09/528,881. 

The magnetic resonance signals from a plurality of magnetic field pulse 
sequences may be used to extract a two-dimensional map of any two parameters 

20 indicative of the fluid in the rock that can be transformed out of the data, such as D, 
Ti, T 2 , viscosity, saturation, etc. The measured data may include thousands of data 
points or, in some cases, tens of thousands of data points or more. In such cases, it 
may be helpful to compress the data before extracting information about the fluid in 
the rock. 
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The two-dimensional density function, /(D,T 2 ), may be extracted from 
Equation 6, 7, or 8, above, using, for example, a least-squares optimization process 
and then plotted to generate a two-dimensional map of diffusion coefficient versus 
relaxation time. Other parameters, such as viscosity, saturation, porosity, etc., may be 
5 derived from D and/or T 2 , and two such parameters mapped against one another to 

generate a two-dimensional map of, for example, T] versus T 2 , viscosity versus 
a saturation, etc. 

Figure 7 shows a full two-dimensional diffusion coefficient-relaxation time 
map, as well as corresponding one-dimensional diffusion and relaxation distributions, 

10 extracted from Equation 7 (A=0.59, B=0.20) using the data from the same partially oil 
saturated sandstone sample discussed above (with reference to Figure 5). The NiCl- 
doped water had about the same relaxation time T 2 as the oil. Thus, in the 
conventional T 2 distribution shown on top (obtained by integrating flJ},T 2 ) over D), 
the oil and water signals overlap completely and only a single peak appears. The two- 

15 dimensional map, as well as the diffusion distribution shown at the right (obtained by 
integrating^^) over T 2 ), clearly show two different diffusing components, one 
with a diffusion coefficient of around 1(T 6 cm 2 /s that can be attributed to the S6 oil 
and the other with a higher diffusion coefficient of around 10" 5 cm 2 /s that can be 
attributed to water. The oil and water saturations correspond to the areas under the 
20 respective peaks, which were found to be 0.78 and 0.22, respectively (in good 
agreement with the results above). 

Thus, according to one aspect of the invention, information such as diffusion 
coefficient and relative saturation of different components of a fluid may be 
determined, even if there is no contrast in relaxation times of the different 
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components. As mentioned before, no assumptions or prior knowledge about the 
relationship between diffusion and relaxation are required to extract this information. 

The two-dimensional map and diffusion distribution in Figure 7 also show 
some contributions having an apparent diffusion coefficient that exceeds that of bulk 
water. These can be attributed to the effects of internal gradients within the rock 
sample which exceed the applied gradient of about 13 Gauss/cm. From Equation 2, it 
can be seen that a large internal gradient may be interpreted as a high diffusion 
coefficient. As mentioned above, the effects of internal gradients can be addressed 
using methods described elsewhere. 

Figure 8 illustrates a plurality of magnetic field pulse sequences that may be 
used in still another embodiment of an NMR measurement according to the invention. 
A first magnetic field pulse sequence 20 has a first portion that includes a first 
stimulated echo sequence. After an initial 90-degree pulse, a second 90-degree pulse 
is applied at time Si followed by a third 90-degree pulse at time Aj,i, stimulating an 
echo at time td T i, all times being relative to the initial 90-degree pulse. A second 
magnetic field pulse sequence 21 is a second stimulated echo-modified CPMG 
sequence in which the time spacings between the initial and the second 90-degree 
pulses, 8 2 , and between the initial and third 90-degree pulses, A2,i, have been modified 
while the time at which a stimulated echo is produced, td,i, is held substantially 
constant. A third and a fourth magnetic field pulse sequences 22 and 23, respectively, 
are also stimulated echo-modified CPMG sequences. In the third magnetic field pulse 
sequence 22, the second 90-degree pulse is generated at a time 5j after the initial 
pulse, as in the first sequence 20, but the time of the stimulated echo, t<j,2, has been 
modified, in this case shortened. The fourth magnetic field pulse sequence 23 
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stimulates an echo at substantially the same time, t d ,2, while the second 90-degree 
pulse is generated at a time 8 2 after the initial pulse, as in the second sequence 21. 

Magnetic resonance signals are detected, as described previously, using the 
second portion of each sequence. By varying die time, td^, at which data processing is 
5 started in addition to varying a time spacing in the first portion in a plurality of 

magnetic field pulse sequences, a time dependent attribute of the fluid in a porous 
* medium, such as a time dependent diffusion coefficient, may be measured. For a 
fluid in a porous medium, such as rock or earth formation, diffusion is restricted. As 
diffusing spins approach the pore wall, the measured diffusion coefficient begins to 
10 deviate from the free fluid diffusion coefficient. As the spins diffuse through several 
pores, the diffusion coefficient asymptotically approaches a tortuosity limit given by 
D 0 /F<p, where D 0 is the free fluid (or bulk) diffusion coefficient, F is the electrical 
formation factor of the porous medium, and <p is the porosity of the medium. Thus, 
from the time dependent behavior of D(t), a direct indication of pore size and/or pore 
15 geometry may be determined. 

Referring to the embodiment illustrated in Figure 8, in the presence of a 
magnetic field gradient, the amplitude of the stimulated echo at time Ui will depend 
on 8j (see Equation 8, above). Thus, measuring how the amplitude of the stimulated 
echo varies with 8j, allows an average diffusion coefficient, D(tda), for all T 2 
20 components at time taj to be extracted, and varying the time tdj for each time spacing 
8j allows a time dependent diffusion coefficient, D(t), to be determined. It should be 
appreciated that other modified CPMG sequences, as described above, may be used in 
time dependent measurements according to the invention, with analogous analyses 
applied. 
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Figure 9 contains graphs showing normalized time dependent diffusion 



coefficient, D(t)/D 0 , versus diffusion length, ^jD 0 t , for two different rock samples. 



The top graph 30 illustrates restricted diffusion in a sandstone sample, while the 
bottom graph 40 illustrates restricted diffusion in a carbonate rock sample. 
5 Photographs of each sample are provided alongside each graph. The NMR 

measurements were conducted with a plurality of stimulated echo-modified CPMG 
sequences at times, Uj, 20 ms, 50 ms, 100 ms, 200 ms, 500 ms, and 1 s. For each 
five values for time spacing 8j were used. The points with error bars, e.g., 32 and 42, 
are the data points generated for each from these measurements. The open circles, 

10 e.g., 34 in the top graph 30 show results from a conventional time dependent diffusion 
coefficient measurement using pulsed field gradients. The dashed lines, 36 and 46, 
indicate the independently determined normalized tortuosity limit, 1/F0 . 

For the sandstone sample, the time dependent diffusion coefficient (as shown 
by points 32) decreases slowly and does not reach the asymptotic tortuosity limit 36 

15 by the longest measured diffusion length (about 50 jam), indicating that the average 
pore size of the sandstone sample is several tens of microns. In contrast, the time 
dependent diffusion coefficient for the carbonate sample (as shown by points 42) 
approaches its tortuosity limit rapidly, indicating an average pore size of a few 
microns. Also, the results of the time dependent diffusion measurements of the 

20 invention on the sandstone sample 32 agree well with the results of conventional 
pulsed field gradient measurements 34. 

Thus, in addition to information such as diffusion and saturation, the magnetic 
field pulse sequences of the invention may be used to extract information about pore 
size and pore geometry of rock without any assumptions or prior knowledge about the 
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rock or the fluid. Conventionally, indications of pore size and pore geometry are 
extracted from relaxation distributions, which depends on phenomenological factors 
such as surface relaxivity, P2, and surface roughness, or from time dependent 
diffusion measurements requiring pulsed field gradients. While the methods of the 
5 invention may be practiced with pulsed field gradients and pulsed field gradients may 
be used to enhance the diffusion effects in a static magnetic field, as discussed 

t previously, pulsed field gradients are not required. The application of pulsed field 

gradients typically requires a large amount of power and additional hardware, such as 
extra RF coils, that in some settings may be undesirable. 

10 The methods of the invention may be practiced in a laboratory setting or in a 

downhole environment, such as with a well logging device. Figure 10 shows an 
apparatus that can be utilized for practicing embodiments of the invention to 
investigate subsurface formations 131 traversed by a borehole 132. A magnetic 
resonance investigating apparatus or logging device 130 is suspended in the borehole 

15 132 on an armored cable 133, the length of which substantially determines the relative 
depth of the device 130. The length of cable 133 is controlled by suitable means at 
the surface such as a drum and winch mechanism. Surface equipment, represented at 
107, can be of conventional type, and can include a processor subsystem that 
communicates with all the downhole equipment. It will be understood that some of 

20. the processing can be performed downhole and that, in some cases, some of the 
processing may be performed at a remote location. Also, while a wireline is 
illustrated, alternative forms of physical support and communicating link can be used, 
for example in a measurement- while-drilling or logging-while-drilling system, in 
practicing the methods of the invention. 
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As described, for example, in U.S. Patent Nos. 5,055,787, 5,055,788, and 
5,153,514, the magnetic resonance logging device 130 can have a face 114 to 
intimately contact the borehole wall. The borehole wall may have a mudcake 116 
thereon. A retractable arm 115 is provided which can be activated to press the body 
5 of the tool 113 through the mudcake against the borehole wall during a logging run, 
with the face 114 pressed against the wall's surface. Although the tool 113 is shown 
as a single body, the tool may alternatively include separate components such as a 
cartridge, sonde or skid, and the tool may be combinable with other logging tools. 

The logging device includes, for example, a permanent magnet or permanent 

10 magnet array 117, which may be made of a samarium-cobalt-magnetic material, and 
one or more RF antennas 118. The investigation region, or sensitivity zone, 
represented generally at 127, is a region in the formation in which the static magnetic 
field is generally uniform, although this is not necessarily required for operation in 
accordance with the invention. Some embodiments of the invention may take 

15 advantage of inherent non- uniformity in the static magnetic field to generate a static 
magnetic field gradient within the investigation region 127. In other embodiments, 
pulsed magnetic field gradients may be used to generate or enhance a magnetic field 
gradient within the investigation region 127. U.S. Patent No. 5,796,252, for example, 
which is incorporated herein by reference, describes various embodiments of an 

20 antenna that can be incorporated into logging devices of the invention and used to 
produce pulse field gradients in the investigation region 127. It will be understood 
that other suitable tool configurations can be utilized for practicing the invention. 

Figure 1 1 shows, in simplified form, circuitry of a type that can be used for 
producing RF pulse sequences and for receiving and processing magnetic resonance 
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signals; it will be appreciated, however, that circuitry having different configurations 
may be used for practicing the invention. 

A downhole processor subsystem is represented at 210. The processor 
subsystem 210 has associated memory, timing, interfaces, and peripherals (not 
5 separately shown), as is known in the art The processor subsystem is conventionally 
coupled with telemetry circuitry 205, for communication with the earth's surface. It 
should be noted that the processor subsystem and associated components may reside 
downhole, uphole, or partially downhole and partially uphole 

An oscillator 220 produces radio frequency (RF) signals at the desired 
10 resonant frequency or frequencies in the investigation region, and the output of the 

oscillator is coupled to a phase shifter 222 and then to a modulator 230, both of which 
are under control of the processor subsystem 210. The phase shifter and modulator 
can be controlled, in a manner known in the art, to produce the desired pulses of RF 
field, for example the 90-degree and 180-degree pulses utilized in embodiments 
15 hereof. As described, for example, in U.S. Patent No. 5,055,788, the oscillator 220 
can be a plurality of oscillators used in a manner that facilitates the generation and 
ultimate detection of the desired signals. The output of modulator 230 is coupled, via 
a power amplifier 235, to the RF antenna 240. A Q-switch 250 can be provided to 
critically dampen the RF antenna system to reduce antenna ringing. 
20 The antenna 240 is also coupled with a receiver section via duplexer 265, the 

output of which is coupled to receiver amplifier 270. The duplexer 265 protects the 
receiver during the transmitting and damping modes. During the receiving mode, the 
duplexer 265 is effectively just a low impedance connection from the antenna to the 
receiver amplifier 270. The output of the receiver amplifier 270 is coupled to a dual 
25 phase-sensitive detector 275, which also receives, as a reference, a signal derived 
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from the oscillator signal. The detected output is coupled to analog-to-digital 
converter 280, the output of which is a digital version of the received nuclear 
magnetic resonance signal. 

Embodiments of the invention may be implemented with logging devices, 
5 such as those described above, without the need for any hardware modifications. 
Figure 12 shows a flow diagram of a routine that can be used in programming a 
processor in implementing embodiments of the invention. The routine may be stored 
on or provided over a computer or machine readable medium, such as read-only 
memory (ROM); random access memory (RAM); magnetic disc or tape; a CD-ROM 
10 or other optical storage media; electrical, optical, acoustical or other forms of 

propagated signals; and the like. The processor may be a downhole processor, an 
uphole processor, or a combination thereof. The processor also may include a remote 
processor that may be used for implementing some of the data interpretation parts of 



magnetic field and, for certain embodiments, a static magnetic field gradient are 
applied to a region of investigation of earth formation surrounding a borehole. 
Theoretically, the earth's magnetic field could be used as the static magnetic field, 
and internal field gradients present within the earth formation could be used as the 

20 static magnetic field gradient. For most practical purposes, however, these are not 
preferred. The static magnetic field and field gradient typically are applied using a 
logging tool having a permanent magnet or an array of permanent magnets, as 
described above. Also, in some embodiments, pulsed magnetic field gradients may be 
applied to the region of investigation from the logging tool. Logging tools can 

25 typically generate static magnetic field gradients in a range of about 10 to 80 



the routine. 



15 



Prior to the beginning of the programmed routine and as shown at 510, a static 
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Gauss/cm, but embodiments of this invention may be implemented with gradients 
outside this range. 

The programmed routine begins at block 515, which represents initializing 
parameters for the magnetic field pulse sequences used in the invention. The 
5 sequence parameters may include, for example, j=l , t^on, td, t long j or 5j and Aj, 

depending on the pulse sequence to be applied. Other sequence parameters may be 
used instead of, or in addition to, the parameters listed (for example, a number (mj) 
may be used to track the number of spin echoes in the first portion, an index, i=l , may 
:: be initialized to track different times etc.), and the sequence parameters may be 
10 initialized all together or at different points in the routine as needed. 

Generating a magnetic field pulse sequence in the region of investigation is 
represented by block 520. The magnetic field pulse sequences used in the invention 
include an initial magnetic field pulse, a first portion and a second portion. Magnetic 
resonance signals from the region of investigation are detected using the second 
15 portion of the magnetic field pulse sequence in block 522 and stored. In some 

embodiments, all magnetic resonance signals from the region of investigation may be 
detected and stored. In other embodiments, magnetic resonance signals preceding a 
time td after the initial pulse of the sequence may not be detected, or signals 
corresponding to time < t<j may be detected but not stored, or stored and later 
20 discarded, or stored and processed in a separate analysis. 

The methods of the invention involve at least two magnetic field pulse 
sequences with differing first portions. A parameter J may be used (set, perhaps, at 
block 515) to indicate a total number of magnetic field pulse sequences to be 
generated. Decision block 525 represents querying whether j = 3. If no, then j is 
25 incremented, as represented in block 526, and sequence parameters, such as t^j, or 5j 



WO01/4281? ^^T/USOO/32589 

^ r 30 

and Aj, or t<j, are varied to modify the first portion of the sequence as represented in 

block 527, before the routine is returned to block 520 where the next sequence is 

generated in the region of investigation. If j=J when queried at decision block 525, 

then the data acquisition for the measurement is complete and the routine continues 

5 on to block 530 where j is re- initialized to j=l. 

In some embodiments of the invention, two or more sequence parameters, for 

example, 8j and t<jj, or t] ongf j and t<jj, are varied sequentially, such that, for each time 

tdj, a plurality of magnetic field pulse sequences are generated in the region of 

investigation using a set of different time spacings, e.g., 5j or t] on gj. Then, time 

10 would be incremented, and a second plurality of magnetic field pulse sequences using 
die same set of time spacings is generated in the region of investigation. To 
implement such embodiments would require initializing an additional parameter, i=l 
(perhaps at block 515), and introducing an additional loop into the routine shown in 
Figure 12 around, for example, blocks 515 and 530. After all the pluralities of 

15 magnetic field pulse sequences have been generated and magnetic resonance signals 
detected and stored, the parameter, i, would be re-initialized to i=l before the routine 
proceeds. 

Blocks 535 and 540 represent, respectively, analyzing the stored magnetic 
resonance signals starting at time U and determining a relative amplitude of the j 111 

20 sequence signal. The relative amplitude of the signal may be quantified as an 

attenuation factor, or the entire data set may be stored for later analysis. Block 545 
represents querying whether all J sequences have been analyzed. If no, parameter j is 
incremented as indicated by block 546 and blocks 535 and 540 repeated until j=J. For 
embodiments of the invention involving an additional set of parameters, such as 

25 times, t^j, an additional loop would be introduced into the routine around, for 
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example, block 535, or around blocks 535 and 545, to analyze the stored signals with 
respect to that parameter. 

When all the relative amplitude data have been calculated, the amplitude data 
is used to extract information about the fluid, as represented by block 550. As 
discussed above, extracting information such as diffusion coefficient, saturation, fluid 
composition, etc. may involve calculating an attenuation factor; or fitting the 
amplitude data to an equation such as Equation 6, 7, or 8 and extracting information 
from the fit; or extracting a full two-dimensional map and identifying different 
components and other information from the map; or combinations of such analyses. 
Additionally, information such as pore size or an indication of pore geometry may be 
extracted from the time dependent behavior of an attribute, such as diffusion 
coefficient, which also can be extracted from the analyzed amplitude data. 

The invention has been described herein with reference to certain examples 
and embodiments. It will, however, be evident that various modifications and 
changes may be made to the embodiments described above without departing from 
the scope of the invention as claimed. 
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We claim: 

1 . A method of extracting information about a fluid in a rock comprising: 

a) generating a sequence of magnetic field pulses in the fluid, the 
sequence comprising an initial magnetic field pulse, a first portion that follows the 
initial magnetic field pulse, and a second portion that follows the first portion; 

b) detecting magnetic resonance signals using the second portion of the 
sequence; 

c) modifying the first portion of the sequence, and repeating steps a) and 

b); and 

d) analyzing magnetic resonance signals from a time relative to the initial 
magnetic field pulse and extracting information about the fluid in the rock. 

2. The method of claim 1, wherein the second portion comprises a series 
of magnetic field pulses separated by a time spacing. 

3. The method of claim 2, wherein the first portion comprises a first 
series of magnetic field pulses separated by a first time spacing. 

4. The method of claim 3, wherein the first time spacing is not less than 
the time spacing of the second portion. 

5. The method of claim 2, wherein the first portion comprises a 
stimulated echo sequence. 
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6. The method of claim 1, wherein the first portion includes at least one 
magnetic field gradient pulse. 

7. The method of claim 1, wherein modifying the first portion comprises 
varying a time spacing between magnetic field pulses. 

8. The method of claim 1, wherein analyzing magnetic resonance signals 
comprises calculating an attenuation factor. 

9. The method of claim 1, wherein analyzing magnetic resonance signals 
comprises determining a relative amplitude of the magnetic resonance signals. 

10. The method of claim 1, wherein extracting information about the fluid 
comprises deteimining a diffusion coefficient. 

1 1 . The method of claim 1 , wherein extracting information about the fluid 
comprises distinguishing between different components of the fluid 

12. The method of claim 1, wherein extracting information about the fluid 
comprises detenniiiing a saturation of the fluid in the rock. 

13. The method of claim 1, further comprising: 
e) repeating step c) a plurality of times. 
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14. The method of claim 13, wherein extracting information about the 
fluid comprises extracting a two dimensional map of a first parameter indicative of 
the fluid in the rock versus a second parameter indicative of the fluid in the rock. 

15. The method of claim 14, wherein the first parameter is a diffusion 
coefficient 

16. The method of claim 14, wherein a second parameter is a relaxation 

time. 



17. The method of claim 13, further comprising: 

f) varying the time relative to the initial magnetic field pulse from which 
the magnetic resonance signals are analyzed and repeating steps a) through e). 

18. The method of claim 17, further comprising repeating step f) a 
plurality of times. 

19. The method of claim 17, wherein extracting information about the 
fluid comprises determining a time dependent diffusion coefficient 

20. The method of claim 17, further comprising determining an indication 
of pore geometry of the rock. 
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21. A metliod of extracting information about a fluid in a rock comprising: 

a) generating a first sequence of magnetic field pulses in the fluid, the 
first sequence comprising a first portion and a second portion; 

b) detecting magnetic resonance signals using the second portion of the 
first sequence; 

c) generating at least one other sequence of magnetic field pulses in the 
fluid, each other sequence comprising a third portion and a fourth portion, the third 

; portion having a different sensitivity to diffusion effects than the first portion of the 
first sequence and the fourth portion having a similar sensitivity to diffusion effects as 
the second portion of the first sequence; 

d) detecting magnetic resonance signals using the fourth portion of each 
other sequence; and 

e) analyzing magnetic resonance signals to separate diffusion effects from 
relaxation effects. 

22. The method of claim 2 1 , wherein the first portion of the first sequence 
is substantially similar to the second portion of the first sequence. 

23. The method of claim 22, wherein the first sequence of magnetic field 
pulses is generated according to: 

90 - [tsho* - 180 - w - echoi]* 
wherein tshoit is a time spacing; echo; is the magnetic resonance spin echo; and n is 
the number of spin echoes. 
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24. The method of claim 21, wherein the first sequence comprises a 
modified CPMG sequence. 

25. The method of claim 24, wherein the first sequence is generated 
according to: 

90 - [ti 0Jlg - 1 80 - tiong - echo k ] ro - [Uh OIi -180- t^on - echor]n- 
wherein t] 0 ji g is a first time spacing; echo* is the magnetic resonance spin echo of 
the first portion; m is the number of spin echoes of the first portion; tshon is a second 
time spacing, the second time spacing being shorter than the first time spacing; echor 
is the i ,th magnetic resonance spin echo of the second portion; and n' is the number of 
spin echoes of the second portion. 

26. The method of claim 24, wherein the first sequence is generated 
according to: 

90 - [5 - 90 - (A-8) - 90 - 5 -echo k ] m - [Wt - 1 80 - t^n - echorV 
wherein 5 and A are time spacings in the first portion; echot is the stimulated 
magnetic resonance spin echo of the first portion; m is the number of stimulated spin 
echoes of the first portion; tshon is a time spacing in the second portion; echo^ is the i' 1 
magnetic resonance spin echo of the second portion; and n' is the number of spin 
echoes of the second portion. 

27. The method of claim 21 , wherein the other sequence comprises a 
modified CPMG sequence. 
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28. The method of claim 27, wherein at least one other sequence is 
generated according to: 

90 - [tj ong j - 1 80 - tjongj ~ echOk j]mj - [tshon - 1 80 - 1^ - echo r j] n - j 
wherein, for the j 01 other sequence, ti OJ1 gj is a third time spacing; echo^j is the k 01 
magnetic resonance spin echo of the third portion; (mj) is the number of spin echoes 
of the third portion; t^on is a fourth time spacing, the fourth time spacing being shorter 
: than the third time spacing; echorj is the V th magnetic resonance spin echo of the 
fourth portion; and (n\j) is the number of spin echoes of the fourth portion. 

29. The method of claim 27, wherein at least one other sequence is 
generated according to: 

90 - [8j - 90 - (Aj-Sj) - 90 - 8j -echokjlmj - [tshon - 180 - t^n - echo r jVj 
wherein, for the j 111 other sequence, Sj and Aj are time spacings in the third portion; 
echo* j is the k* stimulated magnetic resonance spin echo of the third portion; (m,j) is 
the number of stimulated spin echoes of the third portion; tshon is a time spacing in the 
fourth portion; echOj* j is the i ,th magnetic resonance spin echo of the fourth portion; 
and (n\j) is the number of spin echoes of the fourth portion. 

30. The method of claim 21, wherein the first portion and the third portion 
include at least one magnetic field gradient pulse. 

3 1 . The method of claim 2 1 , wherein analyzing magnetic resonance 
signals comprises deter minin g an amplitude of the signal from the other sequence 
relative to an amplitude of the signal from the first sequence. 
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32. The method of claim 21, wherein analyzing magnetic resonance 
signals comprises analyzing spin echoes from the first sequence that correspond in 
time to spin echoes from the other sequence. 

33. The method of claim 21, wherein analyzing magnetic resonance 
signals comprises analyzing a substantially equal number of spin echoes from the first 
sequence and from the other sequence. 

34. The method of claim 21, further comprising extracting a two- 
dimensional map of two parameters indicative of the fluid in the rock. 

35. The method of claim 21, further comprising generating a static 
magnetic field gradient in the rock. 

36. The method of claim 35, further comprising: 

f) moving the rock and the static magnetic field gradient relative to each 



an attribute through at least part of the rock. 

37. The method of claim 21, wherein analyzing magnetic resonance 
signals comprises correlating spin echoes from the first sequence in time with spin 
echoes from the other sequence and analyzing correlated spin echoes at a given time, 
and further comprising: 



other; and 



g) 



repeating steps a) through f) a plurality of times to obtain a profile of 
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f) varying the given time at which correlated spin echoes are analyzed, 

and 

g) repeating steps a) through f) a plurality of times. 

38. The method of claim 37, further comprising determining a time 
dependent attribute of the fluid in the rock. 

39. The method of claim 37, further comprising detenrdning an indication 
of a pore geometry of the rock. 

40. A method of extracting information about a fluid in rock comprising: 

a) generating a first sequence of magnetic field pulses in the fluid, the 
first sequence comprising a first series of magnetic field pulses with a first time 
spacing; 

b) detecting spin echoes using the first series; 

c) generating a second sequence of magnetic field pulses in the fluid, the 
second sequence comprising a second series of magnetic field pulses with a second 
time spacing and a third series of magnetic field pulses with the first time spacing, the 
second spacing being greater than the first time spacing; 

d) detecting spin echoes using the third series; and 

e) extracting information about the fluid in rock using spin echoes 
detected using the first series and spin echoes detected using the third series. 



X:iD: <WO 0142817A1J_> 



40 



7US0Q/32589 



41 . The method of claim 40, wherein extracting information comprises 
analyzing spin echoes detected using the first series that correspond in time to spin 
echoes detected using the third series. 

42. The method of claim 40, wherein extracting information comprises 
analyzing a substantially equal number of spin echoes detected using the first series 
and the third series. 

43. The method of claim 40, wherein extracting information comprises 
determining an amplitude of the spin echoes detected using the third series relative to 
an amplitude of the spin echoes detected using the first series. 

44. The method of claim 40, further comprising: 

f) generating at least one other sequence of magnetic field pulses in the 
fluid, each other sequence comprising a fourth series of magnetic field pulses with a 
time spacing greater than the first time spacing and different from the second time 
spacing and from each other, and a fifth series of magnetic field pulses with the first 
time spacing; and 

g) detecting spin echoes using the fifth series, wherein information about 
the fluid is extracted using the spin echoes detected using the first series, the third 
series, and the fifth series. 
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45. The method of claim 44, wherein extracting information comprises 
analyzing spin echoes detected using the first series that correspond in time to spin 
echoes detected using the fifth series. 

46. The method of claim 44, wherein extracting information comprises 
analyzing a substantially equal number of spin echoes detected using the first series 
and the fifth series. 

47. The method of claim 44, wherein extracting information comprises 
extracting a two-dimensional function of a first parameter and a second parameter 
indicative of the fluid in rock. 

48. The method of claim 47, wherein the first parameter is a diffusion 
coefficient and the second parameter is a relaxation time. 

49. A logging apparatus comprising: 

a logging tool that is moveable through a borehole; and 
a processor that is coupled with the logging tool, the processor being 
programmed with instructions which, when executed by the processor: 
cause the logging tool to: 

i) generate a sequence of magnetic field pulses in a region of 
investigation of earth formation surrounding the borehole, the sequence 
comprising an initial magnetic field pulse, a first portion, and a second 
portion; 
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ii) detect magnetic resonance signals produced from the region of 
investigation using the second portion of the sequence; 

iii) modify the first portion of the sequence and repeat steps (i) and 
(ii); and 

cause the processor to: 

iv) analyze magnetic resonance signals from a time relative to the 
initial magnetic field pulse and extract information about the region of 
investigation. 

50. The logging apparatus of claim 49, wherein the second portion 
comprises a series of magnetic field pulses separated by a time spacing. 

5 1 . The logging apparatus of claim 49, wherein the first portion comprises 
a first series of magnetic field pulses separated by a first time spacing. 

52. The logging apparatus of claim 49, wherein the first portion comprises 
a stimulated echo sequence. 

53. The logging apparatus of claim 49, wherein the instructions, when 
executed by the processor, cause the logging tool to: v) repeat step iii) a plurality 
of times. 



54. The logging apparatus of claim 53, wherein the instructions, when 
executed by the processor, cause the processor to extract a two-dimensional function 
of two parameters that characterize the region of investigation:. 
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55. The logging apparatus of claim 53, wherein the instructions, when 
executed by the processor, cause the processor to vary the time relative to the initial 
magnetic field pulse from which the magnetic resonance signals are analyzed and 
repeat steps i) through v). 

56. A logging apparatus comprising: 

a logging tool that is moveable through a borehole; and 
a processor that is coupled with the logging tool, the processor being 
programmed with instructions which, when executed by the processor: 
cause the logging tool to: 

i) generate a first sequence of magnetic field pulses in a region of 
investigation of earth formation surrounding the borehole, the first sequence 
comprising a first portion and a second portion; 

ii) detect magnetic resonance spin echoes produced from the 
region of investigation using the second portion of the first sequence; 

iii) generate at least one other sequence of magnetic field pulses in 
the region of investigation, each other sequence comprising a third portion and 
a fourth portion, the third portion having a different sensitivity to diffusion 
effects than die first portion of the first sequence and the fourth portion having 
a similar sensitivity to diffusion effects as the second portion of the first 
sequence; 

iv) detect magnetic resonance spin echoes produced from the 
region of investigation using the fourth portion of each other sequence; and 

Cause the processor to: 
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v) analyze detected spin echoes to separate diffusion effects from 
relaxation effects. 



57. The logging apparatus of claim 56, wherein the first portion of the first 
sequence is substantially similar to the second portion of the first sequence. 

58. The logging apparatus of claim 57, wherein the first sequence 
comprises a CPMG sequence. 

59. The logging apparatus of claim 56, wherein the first sequence 
comprises a modified CPMG sequence. 

60. The logging apparatus of claim 56, wherein the other sequence 
comprises a modified CPMG sequence. 

61. The logging apparatus of claim 56, wherein the instructions, when 
executed by the processor, cause the processor to analyze a substantially equal 
number of spin echoes detected using the second portion of the first sequence and 
detected using the fourth portion of the other sequence. 

62. The logging apparatus of claim 56, wherein the instructions, when 
executed by the processor, cause the processor to analyze spin echoes detected using 
the second portion that correspond in time to spin echoes detected using the fourth 
portion. 
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63. The logging apparatus of claim 56, wherein step v) involves analyzing 
spin echoes at a given time, and the instructions, when executed by the processor, 
cause the processor to vary the given time. 

64. A logging apparatus comprising: 

means for generating a sequence of magnetic field pulses in the region of 
investigation, the sequence comprising an initial magnetic field pulse, a first portion 
and a second portion; 

means for detecting magnetic resonance signals using the second portion; 

means for modifying the first portion; and 

means for analyzing magnetic resonance signals from a time relative to the 
initial magnetic field pulse to extract information about the region of investigation. 

65. The logging apparatus of claim 64, further comprising means for 
generating magnetic field gradient pulses in the region of investigation. 

66. A logging apparatus comprising: 

means for generating a first sequence of magnetic field pulses in the region of 
investigation, the first sequence comprising a first portion and a second portion; 

means for detecting magnetic resonance signals from the region of 
investigation using the second portion of the first sequence; 

means for generating at least one other sequence of magnetic field pulses in 
the region of investigation, each other sequence comprising a third portion and a 
fourth portion, the third portion having a different sensitivity to diffusion effects than 
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the first portion of the first sequence and the fourth portion having a similar sensitivity 
to diffusion effects as the second portion of the first sequence; 



investigation using the fourth portion of each other sequence; and 

means for analyzing magnetic resonance signals to separate diffusion effects 
from surface and bulk relaxation effects. 

67. The logging apparatus of claim 66, further comprising means for 
generating magnetic field gradient pulses in the region of investigation. 



means for detecting magnetic resonance signals from the region of 
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NUCLEAR MAGNETIC RESONANCE METHOD AND LOGGING APPARATUS 

This patent application claims priority from U.S. Provisional Application No. 
60/170,121 filed on December 10, 1999 and from U.S. Application Serial No. 
5 09/528,881 filed on March 20, 2000, both of which are herein incorporated by 
reference in their entireties. 

; FIELD OF THE INVENTION 

This invention relates to investigations of rock samples, and more particularly 
10 relates to nuclear magnetic resonance (NMR) methods for determining characteristics 
of subsurface rock, including fluid composition. 

BACKGROUND 

NMR has been a common laboratory technique for over forty years and has 
15 become an important tool in formation evaluation. General background of NMR well 
logging can be found, for example, in U.S. Patent 5,023,551 to Kletnberg et aL, which 
is assigned to the same assignee as the present invention and herein incorporated by 
reference in its entirety. 

NMR relies upon the fact that the nuclei of many chemical elements have 
20 angular momentum ("spin") and a magnetic moment. In an externally applied static 
magnetic field, the spins of nuclei align themselves along the direction of the static 
field. This equilibrium situation can be disturbed by a pulse of an oscillating 
magnetic field (e.g., an RF pulse) that tips the spins away from the static field 
direction. The angle through which the spins are tipped is given by 9=vB]tp/2, where 
25 y is the gyfomagnetic ratio, Bi is the linearly polarized oscillating field strength, and t p 
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is the duration of the pulse. Tipping pulses of ninety and one hundred eighty degrees 
are most common. 

After tipping, two things occur simultaneously. First, the spins precess around 
the direction of the static field at the Larmor frequency, given by coo=yBo, where B 0 is 
5 the strength of the static field and y is the gyromagnetic ratio. For hydrogen nuclei, 
Y/2ti=4258 Hz/Gauss, so, for example, in a static field of 235 Gauss, the hydrogen 
spins would precess at a frequency of 1 MHz. Second, the spins return to the 
equilibrium direction according to a decay time, Ti, which is known as the spin-lattice 
relaxation time. 

10 Also associated with the spin of molecular nuclei is a second relaxation time, 

T2, called the spin-spin relaxation time. At the end of a ninety-degree tipping pulse, 
all the spins are pointed in a common direction perpendicular, or transverse, to the 
static field, and they all precess at the Larmor frequency. However, because of small 
fluctuations in the static field induced by other spins or paramagnetic impurities, the 

15 spins precess at slightly different frequencies, and the transverse magnetization 
dephases with a time constant T 2 . 

A standard technique for measuring T 2 , both in the laboratory and in well 
logging, uses an RF pulse sequence known as the CPMG (Carr-Purcell-Meiboom- 
Gill) sequence. As is well known, after a wait time that precedes each pulse 

20 sequence, a ninety degree pulse causes the spins to start precessing. Then, a one 

hundred eighty degree pulse is applied that keeps the spins in the measurement plane, 
but causes the spins, which are dephasing in the transverse plane, to reverse direction 
and to refocus. By repeatedly reversing the spins using a series of one hundred eighty 
degree pulses, a series of "spin echoes" appear. The train of echoes is measured and 

25 processed to determine the irreversible dephasing, T 2 . 
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In rock formations, such as in a borehole environment, T2 for hydrogen- 



containing fluids (e.g., water, oil, gas) can have significant contributions due to 



surface relaxation, bulk relaxation, and diffusion effects, i.e., 




(1). 



5 



Each of these contributions provides some information about the rock formation 



and/or about the fluid in the rock formation. For example, in a wetting phase, the 
surface relaxation contribution, T2,sujf a ce, dominates the distribution of observed 
distribution of decay times, /(To). Spins relax predominantly due to collisions with 



10 This means that the observed relaxation time is roughly proportional to the pore size, 
i.e., 1 /T 2 ,surface = p2S/V p , where S is the surface area of the pore, V p is the pore volume, 
and p 2 is the surface relaxivity of the rock, a phenomenological parameter that 
indicates how relaxing the surface is. Thus, for a wetting phase, the observed /(T2) 
provides information about pore size distribution. In a nonwetting phase, surface 

15 relaxation becomes negligible and bulk relaxation, which is related to viscosity, 
dominates die observed /(T 2 ). Thus, for a nonwetting phase, the observed J(T 2 ) 
provides information about viscosity. 

In a uniform static magnetic field, each spin will experience the same 
magnetic field strength regardless of its position within the static field, and diffusion 

20 will not contribute to the observed J(T 2 ). In a magnetic field gradient, however, each 
spin will experience different magnetic field strengths as it diffuses through the static 
field. The Larmor frequencies of the diffusing spins become time dependent, and the 
series of one hundred eighty degree pulses cannot refocus the spins completely, 
leading to an additional decay signal. This additional decay signal is proportional to 



the grain surface, with the collision rate being inversely proportional to the pore size. 
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the diffusion coefficient, D, of the fluid and to the square of the gradient strength, g, 
and the square of the echo spacing, t£, i.e., 

■zr^— = jzY 2 8 2 Dtl (2). 

1 2. diffusion i/L 

As the diffusion coefficient provides an indication of fluid type, measurement of the 
5 diffusion effects onyfT 2 ) can be used as the basis for determining the types of fluids in 
a rock formation. 

Certain NMR measurements of diffusion involve changing the echo spacing, 
t E , in a standard CPMG sequence, and thus the amount of diffusion the spins undergo 
between echoes, and then comparing the measured relaxations. Figures 1 A and IB 

10 generally illustrate this approach. Figure 1A shows two CPMG sequences with 
different echo spacings, ti and t 2 , where t 2 is longer than t\. As the echo spacing 
increases, the spins diffuse further between echoes, and the measured relaxation times 
will decrease depending on the diffusion coefficient of the fluid, as given in Equation 
2 above. Figure IB shows the relaxation distributions, /(T 2 ), for an oil and water 

15 determined from the two sets of echoes acquired from the two CPMG sequences 
illustrated in Figure 1A. As seen in Figure IB, the relaxation distribution with the 
longer echo spacing, t 2 , is shifted to lower relaxation times, T 2 , relative to the 
relaxation distribution with the shorter echo spacing, tj. The size of the shift is 
proportional to the size of the diffusion coefficient, as indicated by arrows 1 and 2. 

20 The shift of/[T 2 ) for a fluid with a small diffusion coefficient 1, such as heavy oil, is 
smaller than the shift for a fluid with a larger diffusion coefficient 2, such as water or 
natural gas. 

. While such NMR diffusion measurements can be useful, they suffer from a 
number of drawbacks. For example, for a given acquisition time, the two CPMG 
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sequences will not have the same number of echoes. The CPMG sequence with 
longer echo spacing will have a fewer number of echoes available, so will suffer from 
lower signal to noise and lower data quality in general. Ld addition, relaxation 
distributions for different fluids often overlap, at least partially, making it difficult to 
5 identify shifts of individual relaxation times. In cases where the diffusion coefficients 
for different fluids are small, the shifts may be difficult to distinguish. Finally, these 
methods cannot separate out the contributions due to diffusion effects from the 
surface and bulk relaxation contributions in the observed relaxation distributions. 
Surface relaxation and diffusion have similar effects on the observed relaxation 
10 distributions, so these methods may provide inaccurate information about the fluid 
and about the rock or earth formation under investigation. 

SUMMARY OF INVENTION 

The invention provides in one aspect a method for extracting information 
15 about a fluid that may be contained in rock or within a portion of earth formation 
surrounding a borehole (as used hereinafter, the term "rock" includes earth, earth 
formation, and a portion of earth formation). For example, certain embodiments of 
the invention involve generating a sequence of magnetic field pulses in the fluid. The 
magnetic field pulse sequence includes an initial magnetic field pulse, a first portion 
20 that follows the initial magnetic field pulse, and a second portion that follows the first 
portion. Magnetic resonance signals are detected using the second portion. The first 
portion of the sequence is then modified, and the magnetic field pulse sequence 
generated and magnetic resonance signals detected using the second portion again. 
The magnetic resonance signals are analyzed from a time relative to the initial 
25 magnetic field pulse, and information about the fluid extracted. 
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A second aspect of the invention provides a logging apparatus for 
implementing embodiments of the inventive measurements in a borehole 
environment. One embodiment of such a logging apparatus includes a logging tool 
that is moveable through a borehole and a processor coupled with the logging tool. 
5 The processor is programmed with instructions which, when executed by the 

processor, cause the logging tool to generate a sequence of magnetic field pulses in a 
region of investigation of earth formation surrounding the borehole. The magnetic 
field pulse sequence includes an initial magnetic field pulse, a first portion, and a 
second portion. The processor causes the logging tool to detect magnetic resonance 

10 signals produced from the region of investigation using the second portion of the 

sequence, and then modify the first portion of the sequence and repeat generating the 
sequence and detecting magnetic resonance signals. The programmed instructions 
also cause the processor to analyze magnetic resonance signals from a time relative to 
the initial magnetic field pulse and extract information about the region of 

15 investigation. 

Further details, features and embodiments of the invention will become more 
readily apparent from the detailed description that follows. 



The invention will be described in more detail below in conjunction with the 
20 following Figures, in which: 

Figures 1 A and IB, taken together, illustrate a NMR measurement and T 2 
distributions obtained therefrom according to the prior art; 

Figures 2A and 2B, taken together, illustrate one embodiment of a NMR 
measurement and T 2 distributions obtained dierefrom according to the invention; 
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Figure 3 shows a modified CPMG sequence that can be used in a NMR 
measurement according to embodiments of the invention; 

Figure 4 is a graph of attenuation versus initial echo spacing for three different 
fluids according to certain embodiments of the invention; 
5 Figure 5 is a graph showing relative signal attenuations for three samples as 

initial echo spacing is varied according certain embodiments of the invention; 
p Figure 6 illustrates a plurality of magnetic field pulse sequences that can be 

t . used in one embodiment of a NMR measurement according to the invention; 

Figure 7 shows two-dimensional plot of diffusion coefficient versus relaxation 
10 time according to certain embodiments of the invention; 

Figure 8 illustrates a plurality of magnetic field pulse sequences that can be 
used in another embodiment of a NMR measurement according to the invention; 

Figure 9 shows two graphs illustrating results from two different rock samples 
of a time dependent diffusion coefficient measurement according to certain 
15 embodiments of the invention; 

Figure 10 is a schematic diagram, partially in block form, of one embodiment 
of a well logging apparatus that can be used in implementing methods according to 
the invention in a borehole environment; 

Figure 1 1 is a block diagram of downhole circuitry that can be used with a 
20 well logging apparatus, such as shown in Figure 10; and 

Figure 12 is a flow diagram of a routine that can be used in programming a 
processor or processors in implementing certain embodiments of the invention. 
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DETAILED DESCRIPTION 

A NMR measurement according to the invention involves generating at least 
two different magnetic field pulse sequences in a fluid and detecting magnetic 
5 resonance signals produced by the different magnetic field pulse sequences. 

Information, such as saturation, diffusion coefficient, viscosity, composition, etc., 
about a fluid, e.g., in a rock sample or in a rock formation surrounding a borehole, is 
extracted by analyzing the different magnetic resonance signals. 

Generally speaking, the magnetic field pulse sequences used in the invention 

1 0 can be characterized as having two portions, a first portion followed by a second 

portion. In a NMR measurement according to the invention, the sensitivity of the first 
portion to diffusion effects in the presence of a magnetic field gradient is modified 
while the second portion remains substantially the same. Magnetic resonance signals 
are detected using the second portion and analyzed. Signals detected using the second 

15 portion, which is not changed, will exhibit substantially the same relaxation 

distribution from one magnetic field pulse sequence to the next, except that the 
amplitude of the signals will have been altered as a result of the first portion having 
been modified. By analyzing how the magnetic resonance signals change as the first 
portion is modified, information about the fluid may be extracted. 

20 Stated another way, in the presence of a magnetic field gradient, the first 

portion of the magnetic field pulse sequences used in the invention are sensitive to 
both relaxation and diffusion effects whereas the second portion has substantially the 
same sensitivity to relaxation effects but reduced sensitivity to diffusion effects. 
Diffusion effects during the first portion introduce an extra decay into the signal 

25 detected using the second portion. The signal detected using the second portion is 
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• thus attenuated, or "diffusion edited" in proportion to die diffusion coefficient of the 

fluid (see Equation 2, above). 

These diffusion effects may be detected in the presence of a static magnetic 
field gradient, or with the use of pulsed field gradients as described, for example, in" 
5 U.S. Patent No. 5,796,252 to Klernberg et aL, which is assigned to the same assignee 
as. the present invention and incorporated herein by reference in its entirety. Pulsed 
field gradients introduced into the first portion of the magnetic field pulse sequences 
of the invention also may be used in conjunction with a static magnetic field to 
enhance these diffusion effects. 
10 Figures 2 A and 2B illustrate one embodiment of a NMR measurement 

according to the invention. After a static magnetic field is generated in a fluid in a 
rock, a first magnetic field pulse sequence 3 and then a second magnetic field pulse 
sequence 4 are generated in the fluid. The first sequence 3 in this embodiment 
(shown at the top of Figure 2A) is a standard CPMG sequence in which an initial 90- 
15 degree pulse is followed, after a time tsh ort , by a series of 1 80-degree pulses separated 
in time by about 2t shorl . A magnetic resonance spin echo appears a time t shon after 
each 1 80-degree pulse, producing a series of magnetic spin echoes with a time 
spacing, t E , S bon, approximately equal to 2t shon . This first sequence may be represented 
as: 

20 90 - [When - 1 80 - t sbojl - echo]]* (3) 

where the time separating the echoes, t E>shojt , is equal to about 2t shon ; echo, is the i* 
magnetic resonance spin echo; and n is the number of spin echoes in the sequence. 

A standard CPMG sequence may be characterized as a magnetic field pulse 
sequence according to the invention in which die first portion is substantially identical 

25 to the second portion. A modified CPMG sequence according to the invention may 
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be thought of as a CPMG sequence in which the first portion has been modified. The 
second portion of the modified CPMG sequence is not changed and so generates a 
CPMG-lIke series of magnetic resonance spin echoes with a time spacing 
approximately equal to tE.shon- 
5 The second sequence 4 (shown at the bottom of Figure 2A) is an embodiment 

of a modified CPMG sequence 4 in which the first few echo spacings of the standard 
CPMG sequence are elongated. An initial 90-degree pulse is followed by a first 
portion containing a first series of 1 80-degree pulses that begins a time, t ]ojlgr i, after the 
initial 90-degree pulse and are separated by about 2t] OJlg( i, where t long ,i is greater than 
10 tshoit- A magnetic resonance spin echo appears at a time ti oxl g,] after each 1 80-degree 
pulse, producing a first series of magnetic resonance spin echoes with a time spacing 
tEj approximately equal to 2tj OJlg j. The first portion is followed a second portion 
containing a second series of 1 80-degree pulses separated in time by about 2tshoit. The 
second series of 1 80-degree pulses begins at a time tshort after the last spin echo of the 
15 first portion and refocuses this last spin echo to produce a second series of magnetic 
resonance spin echoes having a time spacing, t£ >S h 0 it, which is equal to about 2tsh<m. 

The embodiment of a modified CPMG sequence 4 shown at the bottom of 
Figure 2A may be represented generally as: 

90 - [tjongj - 1 80 - t Jongj - echoj^j - [^0^ - 180 - tsho.1 - echor j (4) 
20 where, for the j* sequence, t] OJlg j is greater than tshon; echotj is the k 03 magnetic 

resonance spin echo of the first portion; (mj) is the number of spin echoes of the first 
portion; echorj is the i' 1 * 1 magnetic resonance spin echo of the second portion; and 
(n\j) is the number of spin echoes of the second portion. As shown in Figure 2A, the 
first portion of the modified CPMG sequence 4 contains two spin echoes, i.e. y 
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(m,j)=2; it will be appreciated, however, tliat the first portion may have other numbers 
of echoes. 

Another embodiment of a modified CPMG sequence that may be used in 
accordance with the invention has a first portion that contains a stimulated echo 
5 sequence. For example, as shown in Figure 3, the first portion includes a series of 
two 90-degree pulses at times 5j and A] after an initial 90-degree pulse, producing a 
stimulated echo at a time (Ai+Sj) after the initial 90-degree pulse. The second 
portion, which follows the first portion, contains a series of 1 80-degree pulses which 
begins a time t^ on after the last stimulated echo, refocusing the stimulated echo to 
10 produce a series of magnetic resonance spin echoes having a time spacing, t E>sh ort, 

which is equal to about 2^^. The overall sequence may be represented generally as: 
90 - [Sj - 90 - (Aj-Sj) - 90 - Sj -echo kJ ] mJ - [ lshort - 1 80 - tdu* - echor j] VJ (5) 
where, for the j* sequence, echo^j is the k & stimulated echo of the first portion; (m,j) 
is the number of stimulated echoes of the first portion having an echo spacing equal to 
15 about (Aj+Sj); and Sj, Aj, t^, and indices (i\j) and (n'j) are as defined above. 

Other embodiments of a modified CPMG sequence that may be used in 
accordance with the invention are described in U.S. Application Serial No. 
09/528,881 , filed on March 20, 2000, which is incorporated herein by reference. 

The graph in Figure 2B represents the relaxation distributions extracted from 
20 the magnetic resonance signals detected using the second portions of the first 

sequence 3 and the second sequence 4 in the presence of a magnetic field gradient. 
The relaxation distributions, ftj 2 ), for an oil and water, shown in Figure 2B are 
identical except that die amplitude of the second sequence signal 4 is attenuated 
relative to the amplitude of the first sequence signal 3. The amount of the signal 
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attenuation is proportional to the size of the diffusion coefficient, as indicated by 
arrows 5 and 6 Thus, the signal attenuation for a fluid with a small diffusion 
coefficient 5, such as oil, is smaller than the signal attenuation for a fluid with a larger 
diffusion coefficient 6, such as water or natural gas. Measuring the relative change in 
5 amplitude between corresponding spin echoes from the first sequence and from the 
second sequence can yield quantitative information about the fluid in the sample 
under investigation. 

Typically, the echo spacing in the second portion, t E ,shoit 5 is chosen to be as 
short as possible to increase the number, n, of spin echoes that can be generated and 

10 detected within a given acquisition time. This, in general, increases the signal-to- 
noise of the measured signal and, in the presence of a magnetic field gradient, reduces 
the sensitivity of the second portion to diffusion effects. Times for tshon on the order 
of about 0. 1 milliseconds (100 fas), leading to echo spacings, tE,shoit, on the order of 
about 0.2 ms (200 jas), are currently used in well logging measurements, although the 

15 measurements of the invention may be made using other times for tshon and t£ >s hoit- 
To properly correlate the spin echoes from the first sequence with the spin 
echoes from the second sequence, the data processing for both sequences starts at the 
same time, t<j, relative to the initial 90-degree pulse. In other words, spin echoes that 
precede t<i are not used in processing, and only spin echoes starting from U are 

20 analyzed. By beginning the data processing at the same time for each sequence, the 
magnetic resonance signals from each sequence will reflect substantially similar 
relaxation effects. 

As shown in Figure 2A, t<j corresponds to the last echo of the first portion of 
. the modified CPMG sequence. For example, with tsh6n=0. 1 ms and ti ong j=4 ms, data 
25 processing would start at t s =16 ms, and if the two spin echoes of the first portion of 
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the second sequence are disregarded, then the first 80 spin echoes of the first 
sequence, i.e., echo; from i=l to (mj)x(ti 0 ngj/t dl0rt ), would be disregarded. Data 
processing according to the certain embodiments of the invention, however, may 
include the last echo of the first portion, i.e., the echo from the first portion that is 
5 re focused by the magnetic pulses of the second portion. This means that the spin 
echoes that occur at t=t<j would be included in the data processing. 

Both sequences also may be truncated at the end during processing so a 
substantially equal number of spin echoes are analyzed and/or echoes with poor 
signal-to-noise are disregarded. Typically, both sequences still will contain a large 
10 number of echoes, particularly as compared to prior art techniques in which the echo 
spacing of a standard CPMG is increased. Thus, the methods of the invention will 
generally allow more precise measurements and provide better signal-to-noise as 
compared to prior art techniques. 

Though the spin echoes preceding td typically are not used in the above 
15 analysis, they may be collected and analyzed to extract other information, such as 
total porosity, permeability,. etc. of a fluid-containing rock sample, as described, for 
example, in U.S. Patent Nos. 5,023,551, 5,363,041, and 5,389,877. Such analyses 
may be done in conjunction with the methods of this invention. 

A diffusion edited signal generated according to (4), for times > ta, may be 
20 represented as: 

M (hj,0 = jjdDdT 1 f(D,T 2 ) & xp(-^)e W (-^- r 2 8 2 Dtl J ) (6), 

where t<j conresponds to the time (relative to the initial pulse) of the last echo of the 
first portion^D/IV) is the two-dimensional diffusion-T 2 probability density function, 
(m j) is the number of echoes of the first portion, y is. the gyromagnetic ratio, g is the 
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gradient field strength, D is the diffusion coefficient, and t E j is the echo spacing in 
the first portion of the modified sequence. For modified sequences having more than 
one echo in the first portion, Le. t (m,j)>l, multi exponential diffusion decays may be 
observed over a broad detection bandwidth. Such multiexponential decay can be 
5 modeled more accurately by replacing the single exponential diffusion attenuation 
term in Equation 6 with multiple terms. For example, for (m,j)=2, two exponential 
terms, one from the direct echo and the other from the stimulated echo, would replace 
the single exponential term in Equation 6, yielding: 

M(/ £J J) = J/JZ^ (7), 

10 where A and B are parameters that depend only on the detection bandwidth 

Similarly, a diffusion edited signal generated according to (5), for times > t<j, 
may be represented as: 

where U is the time (relative to the initial pulse) of the stimulated echo,/(D(td),T2) is 
15 the two-dimensional diffusion-T 2 probability density function at time td, and D(t<0 
represents the diffusion coefficient at time t<j. 

The amplitude of the detected echoes can be measured using any of various 
signal processing techniques known in the art and then, according to certain 
embodiments of the invention, fit to one of Equation 6, 7, 8, or other diffusion edited 
20 signal equations, depending on the magnetic pulse sequence used, in order to extract 
diffusion coefficient and other information about the fluid. 

As previously suggested, measuring the relative change in amplitude between 
corresponding spin echoes from the first and second magnetic field pulse sequences 
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cao provide quantitative information about the fluids in the sample under 
investigation. In particular, a comparison of corresponding spin echoes can separate 
different contributions to the T 2 relaxation distribution, and so yield more accurate 
information about the fluid. For example, as discussed with regards to Figure 2A, in 
5 the presence of a magnetic field gradient, the first portion of the modified CPMG 
sequences of the invention are more sensitive to diffusion effects than the second 
2 portion, leading to the signal from the second portion to be attenuated, or diffusion- 
edited, compared to the standard CPMG signal. However, the signal from the second 
portion exhibits the same surface and bulk relaxation effects as, and is otherwise 

10 identical to, the signal from the standard CPMG sequence. Thus, taking a ratio of the 
amplitudes of the corresponding T 2 distributions of the first and second sequences 
separates diffusion effects from both surface and bulk relaxation effects. 

For multiphase fluids, linear mixing laws govern the observed signal 
attenuation. Thus, components of the fluid having different diffusion coefficients 

15 may be differentiated by appropriately selecting and varying a time spacing in the first 
portion of the modified sequences (e.g., t£,i; Sj or Aj). Figure 4 contains a graph that 
may be useful in selecting a first portion time spacing. The graph shows approximate 
attenuation versus initial echo spacing, t E j, for oil (having a diffusion coefficient, D 0 a, 
of 2.3 x 10' 6 cm 2 /s), water (D walei =2.3 x 10" 5 cm 2 /s) and gas (D gas =1.2 x 10" 3 cm 2 /s) 

20 measured in a magnetic field gradient of about 25 Gauss/cm with a modified sequence 
generated according to (4), having two echoes in the first portion. Over a range of 
echo spacings greater than approximately 0.8 ms, the relative attenuation of the gas 
signal is much larger than for either water or oil, indicating that this range may be 
useful in differentiating gas. from water or oil. In general, an initial echo spacing in a 
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10 



15 



range of between about 0.5 to about 20 ms may be useful for hydrocarbon typing 
applications of the invention. 

An initial echo spacing to differentiate between two diffusing components 
may be approximated as: 



where D + is the diffusion coefficient of the higher diffusing component, D. is the 
diffusion coefficient of the lower diffusing component, y is the gyromagnetic ratio and 
g is the magnetic field gradient strength. Equation 9 may provide a useful starting 
point in selecting an initial time spacing, if the different diffusing components (or a 
range of diffusion coefficients) are known or can be approximated a priori. 

In some embodiments of the invention, the relative attenuation of the echo 
signal may be quantified in terms of an attenuation factor, a;. The fluid composition 
then may be determined by taking a linear combination of the attenuation factors for 
the individual components. For example, the attenuation factor may be calculated as 
a ratio of the sums of all the measured echoes of a modified CPMG sequence to the 
measured echoes of a standard CPMG sequence, or as a ratio of initial amplitudes of a 
modified CPMG sequence to a standard CPMG sequence. In a gradient field of about 
25 Gauss/cm and using a modified sequence generated according to (4) with two 
echoes in the first portion (t E j = 8 ms, t E ,short = 0.2 ms), an attenuation factor for bulk 
water, a w , was calculated to be about 0.32 and for 6 cp oil, aou was calculated to be 
about 0.85. For a fluid having water and oil components, the water saturation, S w , i.e., 
the proportion of fluid that is water, can be determined from the measured attenuation 
factor of the fluid, ameas, u sing the following relationship: 



t 




(9), 



EJ ~ 



Y 2 g 2 (D+-D_) 
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^ _ ~" ^mgax (10). 

According to other embodiments, the fluid composition can be determined by 
fitting the measured attenuation curve versus a first portion time spacing as a 
superposition of component curves. Figure 5 shows the relative initial attenuations of 
5 diffusion edited signals as an initial echo spacing is changed for a NiCl-doped water 
sample (labeled "water"), a S6 oil sample (labeled "oil"), and a partially oD saturated 
sandstone sample. The diffusion edited signals were generated in a uniform gradient 
field of about 13 Gauss/cm with a modified CPMG sequence according to (4) and 
having two echoes in the first portion. For the water and oil samples, the plotted 

10 points, representing the amplitude of a modified CPMG signal relative to a standard 
CPMG signal, were fit to a curve using Equation 7 above, with A=0.59 and B=0.20. 
From the fit, the diffusion coefficients were extracted and found to be Dw=£ xlO 5 
cm 2 /s for the water sample and Dou=1.35xl0" 6 cm 2 /s for the oil sample, "or he 
partially oil saturated sandstone sample, the relative attenuation was modeled as a 

15 superposition of two curves using Equation 7 and three parameters: S w , the* ater 

saturation; D w ,eff, the effective diffusion coefficient for water; and D oi i >e ff. the : ffective 
diffusion coefficient for oil. From the fit, it was determined that: S w =0. Z\ , 
D Wt efF=1.9xl0" 5 cm 2 /s, and D O ii rC ff=13xl0" 6 - cm 2 /s. 

The water saturation, S w , determined in this manner agreed with independent 

20 gravimetrical measurements, and the effective oil diffusion coefficient D ... c matched 
the bulk diffusion coefficient D c y of the pure fluid, as expected. The effective- water 
diffusion coefficient D w ,eff, however, was smaller than the bulk diffusioi. cot. ficient, 
D w , of water. This difference in water diffusion coefficients is attributed to restricted 



)OClD: <WO 0142817A1JA> 



. WO 01/42817 




7US00/32589 



diffusion of the water in the rock, as water in sandstone is known to mainly occupy 
the small pores of the rock. 

The effects of pore geometry, in terms of restricted diffusion and internal 
gradients, may be addressed using the methods described in U.S. Application Serial 
5 No. 09/528,881. The effects of restricted diffusion and internal gradients may also be 
addressed in terms of the measured attenuation factor. A hypothetical water 
saturation for water that does not contribute to the observed signal attenuation, 

calculated as, e.g., S w regtr = a,neas - ° w , 

and a ratio of the measured attenuation factor to that of bulk water, i.e., Qfn<?a5 , may be 

10 used to correct for the effects of restricted diffusion and internal gradients, 

respectively, in calculating, for example, water/oil saturation in a rock sample. 

The diffusion editing embodiments of the invention thus can extract saturation 
and diffusion coefficient information about a fluid directly without a priori knowledge 
of, and without having to make any assumptions about, relaxation in the fluid. As the 

15 measurements of these embodiments take place in the presence of a magnetic field 
gradient, the spin echoes are generated from a thin slice across the sample. By 
moving the sample and the gradient relative to one another along the gradient 
direction, an attribute profile of the fluid, for example, a saturation profile of a fluid in 
a rock sample, can be obtained. 

20 In other embodiments of the invention, a plurality of magnetic field pulse 

sequences, each with a different first portion (e.g., different time spacing, t^j; or 5j 
and/or Aj), are used in making the NMR measurements. The diffusion edited 
magnetic resonance signal depends on a time spacing in the first portion (e.g., t E j; 5j, 
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Aj), so fluids with different diffusion coefficients will diffuse different amounts during 
different time spacings (see Equations 6, 7, or 8, above). Inverting data using only 
two different time spacings (as described in the embodiments above) yields an 
average diffusion coefficient for every relaxation time and so may not resolve 
5 different diffusing components that have overlapping relaxation times. Measuring 
magnetic resonance signals using more than two different time spacings allows 
diffusion coefficient and relaxation time to be extracted separately and so can help 
resolve different fluid components in terms of diffusion and relaxation, as well as 
other parameters that may be derived from D, T 2 , or a combination of D and T 2 , such 

10 as Ti, viscosity, saturation, porosity, etc.. 

Figure 6 illustrates a plurality of magnetic field pulse sequences that may be 
used in one embodiment of a NMR measurement according to the invention. A first 
magnetic field pulse sequence 10 is a standard CPMG sequence having a first time 
spacing, t^on, and a first echo spacing, t EtS hoit- A second magnetic field pulse sequence 

15 11 is generated by modifying, in this case elongating, a time spacing between the 
initial 90-degree pulse and the first 180-degree pulse, ti ong ,i, leading to an elongated 
echo spacing, t E ,i, in the first portion. A third magnetic field pulse sequence 12 is 
generated by modifying a time spacing of the first portion again, t] ong> 2- A fourth 
magnetic, field pulse sequence 13 is again generated by modifying a time spacing of 

20 the first portion, ti on g,3, and also the number of spin echoes in the first portion. 

This embodiment of a NMR measurement involves generating each of the 
magnetic field pulse sequences shown in Figure 6 in a fluid in a rock. Magnetic 
resonance signals from each sequence are detected using the second portion of each 
sequence, which remains substantially the same from sequence to sequence. As 

25 discussed above, to properly correlate the spin echoes detected from each sequence, 
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the data processing for each sequence starts at a time u relative to the initial 90-degree 
pulse. Figure 6 shows t<j corresponding to the time of the last echo of the second 
sequence 11, which corresponds in this illustration with the time of the last echo of 
the fourth sequence 13. Some spin echoes from the second portion of the third 
5 sequence 12, which arise prior to time td, are discarded in the data processing 

according to this embodiment The time u need not necessarily correspond with the 
time of the last echo of the longest first portion, as shown in Figure 6; however, t<j 
typically is at least as long as the longest first portion. 

Although Figure 6 shows four magnetic field pulse sequences, other numbers 

10 of magnetic field pulse sequences may be used in NMR measurements according to 
the invention, with the use of more sequences resulting in higher resolution in the 
diffusion (or relaxation, viscosity, saturation, etc.) distribution. Additionally, the 
plurality of magnetic field pulse sequences need not include one type of modified 
CPMG sequence, as shown in Figure 6. Measurements according to methods of the 

15 invention may be made using a plurality of magnetic field pulse sequences that 

includes a combination of modified CPMG sequences generated according to (4), (5) 
or any of the other sequences described in U.S. Application Serial No. 09/528,881. 

The magnetic resonance signals from a plurality of magnetic field pulse 
sequences may be used to extract a two-dimensional map of any two parameters 

20 indicative of the fluid in the rock that can be transformed out of the data, such as D, 
Ti, To, viscosity, saturation, etc. The measured data may include thousands of data 
points or, in some cases, tens of thousands of data points or more. In such cases, it 
may be helpful to compress the data before extracting information about the fluid in 
the rock. 
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The two-dimensional density function, y(D,T 2 ), may be extracted from 
Equation 6, 7, or 8, above, using, for example, a least-squares optimization process 
and then plotted to generate a two-dimensional map of diffusion coefficient versus 
relaxation time. Other parameters, such as viscosity, saturation, porosity, etc., may be 
5 derived from D and/or T 2 , and two such parameters mapped against one another to 
generate a two-dimensional map of, for example, Ti versus To, viscosity versus 
saturation, etc. 

Figure 7 shows a full two-dimensional diffusion coefficient-relaxation time 
map, as well as corresponding one-dimensional diffusion and relaxation distributions, 

10 extracted from Equation 7 (A=0.59, B=0.20) using the data from the same partially oil 
saturated sandstone sample discussed above (with reference to Figure 5). The NiCl- 
doped water had about the same relaxation time T 2 as the oil. Thus, in the 
conventional T 2 distribution shown on top (obtained by integrating y(D,T 2 ) over D), 
the oil and water signals overlap completely and only a single peak appears. The two- 

15 dimensional map, as well as the diffusion distribution shown at the right (obtained by 
integrating y(D,T 2 ) over T 2 ), clearly show two different diffusing components, one 
with a diffusion coefficient of around 1CT 6 cm 2 /s that can be attributed to the S6 oil 
and the other with a higher diffusion coefficient of around 10" 5 cm 2 /s that can be 
attributed to water. The oil and water saturations conrespond to the areas under the 

20 respective peaks, which were found to be 0.78 and 0.22, respectively (in good 
agreement with the results above). 

Thus, according to one aspect of the invention, information such as diffusion 
coefficient and relative saturation of different components of a fluid may be 
determined, even if there is no contrast in relaxation times of the different 
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components. As mentioned before, no assumptions or prior knowledge about the 
relationship between diffusion and relaxation are required to extract this information. 

The two-dimensional map and diffusion distribution in Figure 7 also show 
some contributions having an apparent diffusion coefficient that exceeds that of bulk 
5 water. These can be attributed to the effects of internal gradients within the rock 

sample which exceed the applied gradient of about 13 Gauss/cm. From Equation 2, it 
can be seen that a large internal gradient may be interpreted as a high diffusion 
coefficient. As mentioned above, the effects of internal gradients can be addressed 
using methods described elsewhere. 
10 Figure 8 illustrates a plurality of magnetic field pulse sequences that may be 

used in still another embodiment of an NMR measurement according to the invention. 
A first magnetic field pulse sequence 20 has a first portion that includes a first 
stimulated echo sequence. After an initial 90-degree pulse, a second 90-degree pulse 
is applied at time Si followed by a third 90-degree pulse at time A u , stimulating an 
15 echo at time td.i, all times being relative to the initial 90-degree pulse. A second 
magnetic field pulse sequence 21 is a second stimulated echo-modified CPMG 
sequence in which the time spacings between the initial and the second 90-degree 
pulses, 82, and between the initial and third 90-degree pulses, A 2 ,i, have been modified 
while the time at which a stimulated echo is produced, is held substantially 
20 constant. A third and a fourth magnetic field pulse sequences 22 and 23, respectively, 
are also stimulated echo-modified CPMG sequences. In the third magnetic field pulse 
sequence 22, the second 90-degree pulse is generated at a time Si after the initial 
pulse, as in the first sequence 20, but the time of the stimulated echo, has been 
modified, in this case shortened. The fourth magnetic field pulse sequence 23 

■ : t \ 
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stimulates an echo at substantially the same time, td, 2 , while the second 90-degree 
pulse is generated at a time 5 2 after the initial pulse, as in the second sequence 21. 

Magnetic resonance signals are detected, as described previously, using the 
second portion of each sequence. By varying the time, t^, at which data processing is 
5 started in addition to varying a time spacing in the first portion in a plurality of 
magnetic field pulse sequences, a time dependent attribute of the fluid in a porous 
medium, such as a time dependent diffusion coefficient, may be measured. For a 
fluid in a porous medium, such as rock or earth formation, diffusion is restricted. As 
diffusing spins approach the pore wall, the measured diffusion coefficient begins to 

10 deviate from the free fluid diffusion coefficient. As the spins diffuse through several 
pores, the diffusion coefficient asymptotically approaches a tortuosity limit given by 
D 0 /F<p , where D 0 is the free fluid (or bulk) diffusion coefficient, F is the electrical 
formation factor of the porous medium, and <p is the porosity of the medium. Thus, 
from the time dependent behavior of D(t), a direct indication of pore size and/or pore 

15 geometry may be determined. 

Referring to the embodiment illustrated in Figure 8, in the presence of a 
magnetic field gradient, the amplitude of the stimulated echo at time t^j will depend 
on Sj (see Equation 8, above). Thus, measuring how the amplitude of the stimulated 
echo varies with 5j, allows an average diffusion coefficient, D(tda), for all T 2 

20 components at time tdj to be extracted, and varying the time t^ for each time spacing 
5j allows a time dependent diffusion coefficient, D(t), to be determined. It should be 
appreciated that other modified CPMG sequences, as described above, may be used in 
time dependent measurements according to the invention, with analogous analyses 
applied. . 
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Figure 9 contains graphs showing normalized time dependent diffusion 



coefficient, D(t)/D Q , versus diffusion lengdi, ^/Iv > f° r two different rock samples. 
The top graph 30 illustrates restricted diffusion in a sandstone sample, while the 



5 Photographs of each sample are provided alongside each graph. The NMR 

measurements were conducted with a plurality of stimulated echo-modified CPMG 
sequences at times, taj, 20 ms, 50 ms, 100 ms, 200 ms, 500 ms, and 1 s. For each t^, 
five values for time spacing 6j were used. The points with error bars, e.g., 32 and 42, 
are the data points generated for each tdj from these measurements. The open circles, 
10 e.g., 34 in the top graph 30 show results from a conventional time dependent diffusion 
coefficient measurement using pulsed field gradients. The dashed lines, 36 and 46, 
indicate the independently determined normalized tortuosity limit, 1/ F(p . 

For the sandstone sample, the time dependent diffusion coefficient (as shown 
by points 32) decreases slowly and does not reach the asymptotic tortuosity limit 36 

15 by the longest measured diffusion length (about 50 Jim), indicating that the average 
pore size of the sandstone sample is several tens of microns. In contrast, the time 
dependent diffusion coefficient for the carbonate sample (as shown by points 42) 
approaches its tortuosity limit rapidly, indicating an average pore size of a few 
microns. Also, the results of the time dependent diffusion measurements of the 

20 invention on the sandstone sample 32 agree well with the results of conventional 
pulsed field gradient measurements 34. 

Thus, in addition to information such as diffusion and saturation, the magnetic 
field pulse sequences of the invention may be used to extract information about pore 
size and pore geometry of rock without any assumptions or prior knowledge about the 



bottom graph 40 illustrates restricted diffusion in a carbonate rock sample. 
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rock or the fluid. Conventionally, indications of pore size and pore geometry are 
extracted from relaxation distributions, which depends on phenomenological factors 
such as surface relaxivity, p 2 , and surface roughness, or from time dependent 
diffusion measurements requiring pulsed field gradients. While the methods of the 
5 invention may be practiced with pulsed field gradients and pulsed field gradients may 
be used to enhance the diffusion effects in a static magnetic field, as discussed 
previously, pulsed field gradients are not required. The application of pulsed field 
gradients typically requires a large amount of power and additional hardware, such as 
extra RP coils, that in some settings may be undesirable. 

10 The methods of the invention may be practiced in a laboratory setting or in a 

downhole environment, such as with a well logging device. Figure 10 shows an 
apparatus that can be utilized for practicing embodiments of the invention to 
investigate subsurface formations 131 traversed by a borehole 132. A magnetic 
resonance investigating apparatus or logging device 130 is suspended in the borehole 

15 132 on an armored cable 133, the length of which substantially determines the relative 
depth of the device 130. The length of cable 133 is controlled by suitable means at 
the surface such as a drum and winch mechanism. Surface equipment, represented at 
107, can be of conventional type, and can include a processor subsystem that 
communicates with all the downhole equipment. It will be understood that some of 

20 the processing can be performed downhole and that, in some cases, some of the 
processing may be performed at a remote location. Also, while a wireline is 
illustrated, alternative forms of physical support and communicating link can be used, 
for example in a measurement- while-drilling or logging- while-drilling system, in 
practicing the methods of the invention. 
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As described, for example, in U.S. Patent Nos. 5,055,787, 5,055,788, and 
5,153,5 14, the magnetic resonance logging device 130 can have a face 114 to 
intimately contact the borehole wall. The borehole wall may have a mudcake 116 
thereon. A retractable arm 115 is provided which can be activated to press the body 
5 of the tool 113 through the mudcake against the borehole wall during a logging run, 
with the face 114 pressed against the wall's surface. Although the tool 113 is shown 
as a single body, the tool may alternatively include separate components such as a 
cartridge, sonde or skid, and the tool may be combinable with other logging tools. 

The logging device includes, for example, a permanent magnet or permanent 

10 magnet array 117, which may be made of a samarium-cobalt-magnetic material, and 
one or more RF antennas 118. The investigation region, or sensitivity zone, 
represented generally at 127, is a region in the formation in which the static magnetic 
field is generally uniform, although this is not necessarily required for operation in 
accordance with the invention. Some embodiments of the invention may take 

15 advantage of inherent non-uniformity in the static magnetic field to generate a static 
magnetic field gradient within the investigation region 127. In other embodiments, 
pulsed magnetic field gradients may be used to generate or enhance a magnetic field 
gradient within the investigation region 127. U.S. Patent No. 5,796,252, for example, 
which is incorporated herein by reference, describes various embodiments of an 

20 antenna that can be incorporated into logging devices of the invention and used to 
produce pulse field gradients in the investigation region 127. It will be understood 
that other suitable tool configurations can be utilized for practicing the invention. 

Figure 1 1 shows, in simplified form, circuitry of a type that can be used for 
producing RF pulse sequences and for receiving and processing magnetic resonance 
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20 



signals; it will be appreciated, however, that circuitry having different configurations 
may be used for practicing the invention. 



subsystem 210 has associated memory, timing, interfaces, and peripherals (not 
separately shown), as is known in the art. The processor subsystem is conventionally 
coupled with telemetry circuitry 205, for communication with the earth's surface. It 
should be noted that the processor subsystem and associated components may reside 
downhole, uphole, or partially downhole and partially uphole 

An oscillator 220 produces radio frequency (RF) signals at the desired 
resonant frequency or frequencies in the investigation region, and the output of the 
oscillator is coupled to a phase shifter 222 and then to a modulator 230, both of which 
are under control of the processor subsystem 210. The phase shifter and modulator 
can be controlled, in a manner known in the art, to produce the desired pulses of RF 
field, for example the 90-degree and 180-degree pulses utilized in embodiments 
hereof. As described, for example, in U.S. Patent No. 5,055,788, the oscillator 220 
can be a plurality of oscillators used in a manner that facilitates the generation and 
ultimate detection of the desired signals. The output of modulator 230 is coupled, via 
a power amplifier 235, to the RF antenna 240. A Q-switch 250 can be provided to 
critically dampen the RF antenna system to reduce antenna ringing. 

The antenna 240 is also coupled with a receiver section via duplexer 265, the 
output of which is coupled to receiver amplifier 270. The duplexer 265 protects the 
receiver during the transmitting and damping modes. During the receiving mode, the 
duplexer 265 is effectively just a low impedance connection from the antenna to the 
receiver amplifier 270. The output of the receiver amplifier 270 is coupled to a dual 
phase-sensitive detector 275, which also receives, as a reference, a signal derived 



A downhole processor subsystem is represented at 210. The processor 
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from the oscillator signal. The detected output is coupled to analog-to-digital 
converter 280, the output of which is a digital version of the received nuclear 
magnetic resonance signal. 

Embodiments of the invention may be implemented with logging devices, 
5 such as those described above, without the need for any hardware modifications. 
Figure 12 shows a flow diagram of a routine that can be used in programming a 
processor in implementing embodiments of the invention. The routine may be stored 
on or provided over a computer or machine readable medium, such as read-only 
memory (ROM); random access memory (RAM); magnetic disc or tape; a CD-ROM 
10 or other optical storage media; electrical, optical, acoustical or other forms of 

propagated signals; and the like. The processor may be a downhole processor, an 
uphole processor, or a combination thereof. The processor also may include a remote 
processor that may be used for implementing some of the data interpretation parts of 
the routine. 

15 Prior to the beginning of the programmed routine and as shown at 510, a static 

magnetic field and, for certain embodiments, a static magnetic field gradient are 
applied to a region of investigation of earth formation surrounding a borehole. 
Theoretically, the earth's magnetic field could be used as the static magnetic field, 
and internal field gradients present within the earth formation could be used as the 

20 static magnetic field gradient. For most practical purposes, however, these are not 
preferred. The static magnetic field and field gradient typically are applied using a 
logging tool having a permanent magnet or an array of permanent magnets, as 
described above. Also, in some embodiments, pulsed magnetic field gradients may be 
applied to the region of investigation from the logging tool. Logging tools can 

25 typically generate static magnetic field gradients in a range of about 10 to 80 



vJSOOCID: <WO 0142817A1 JA> 



.WO 01/42817 \ ^P/US00/32589 

29 

Gauss/cm, but embodiments of this invention may be implemented with gradients 
outside this range. 

The programmed routine begins at block 515, which represents initializing 
parameters for the magnetic field pulse sequences used in the invention. The 
5. sequence parameters may include, for example, j=l , LW, td, t Jojlg j or 5j and Aj, 

depending on the pulse sequence to be applied. Other sequence parameters may be 
used instead of, or in addition to, the parameters listed (for example, a number (m,j) 
may be used to track the number of spin echoes in the first portion, an index, i=l, may 
be initialized to track different times taj, etc.), and the sequence parameters may be 
10 initialized all together or at different points in the routine as needed. 

Generating a magnetic field pulse sequence in the region of investigation is 
represented by block 520. The magnetic field pulse sequences used in the invention 
include an initial magnetic field pulse, a first portion and a second portion. Magnetic 
resonance signals from the region of investigation are detected using the second 
15 portion of the magnetic field pulse sequence in block 522 and stored. In some 

embodiments, all magnetic resonance signals from the region of investigation may be 
detected and stored. In other embodiments, magnetic resonance signals preceding a 
time u after the initial pulse of the sequence may not be detected, or signals 
corresponding to time < u may be detected but not stored, or stored and later 
20 discarded, or stored and processed in a separate analysis. 

The methods of the invention involve at least two magnetic field pulse 
sequences with differing first portions. A parameter J may be used (set, perhaps, at 
block 515) to indicate a total number of magnetic field pulse sequences to be 
generated. Decision block 525 represents querying whether j = J. If no, thenj is 
25 incremented, as represented in block 526, and sequence parameters, such as t] on gj, or 5j 
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and Aj, or U, are varied to modify the first portion of the sequence as represented in 
block 527, before the routine is returned to block 520 where the next sequence is 
generated in the region of investigation. If j=J when queried at decision block 525, 
then the data acquisition for the measurement is complete and the routine continues 
5 on to block 530 where j is re- initialized to j=l. 

In some embodiments of the invention, two or more sequence parameters, for 
example, 8j and taa, or tiang.j and t^, are varied sequentially, such that, for each time 
t<jj, a plurality of magnetic field pulse sequences are generated in the region of 
investigation using a set of different time spacings, e.g., 8j or t] 0n gj- Then, time tdj 

10 would be incremented, and a second plurality of magnetic field pulse sequences using 
the same set of time spacings is generated in the region of investigation. To 
implement such embodiments would require initializing an additional parameter, i=l 
(perhaps at block 515), and introducing an additional loop into the routine shown in 
Figure 12 around, for example, blocks 515 and 530. After all the pluralities of 

15 magnetic field pulse sequences have been generated and magnetic resonance signals 
detected and stored, the parameter, i, would be re-initialized to i=l before the routine 
proceeds. 

Blocks 535 and 540 represent, respectively, analyzing the stored magnetic 
resonance signals starting at time U and deteniiining a relative amplitude of the j 111 

20 sequence signal. The relative amplitude of the signal may be quantified as an 

attenuation factor, or the entire data set may be stored for later analysis: Block 545 
represents querying whether all J sequences have been analyzed. If no, parameter j is 
incremented as indicated by block 546 and blocks 535 and 540 repeated until j=J. For 
embodiments of the invention involving an additional set of parameters, such as 

25 times, t^, an additional loop would be introduced into the routine around, for 
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example, block 535, or around blocks 535 and 545, to analyze the stored signals with 

respect to that parameter. 

When all the relative amplitude data have been calculated, the amplitude data 

is used to extract information about die fluid, as represented by block 550. As 
5 discussed above, extracting information such as diffusion coefficient, saturation, fluid 

composition, etc. may involve calculating an attenuation factor; or fitting the 
t amplitude data to an equation such as Equation 6, 7, or 8 and extracting information 

from the fit; or extracting a full two-dimensional map and identifying different 

components and other information from the map; or combinations of such analyses. 
10 Additionally, information such as pore size or an indication of pore geometry may be 

extracted from the time dependent behavior of an attribute, such as diffusion 

coefficient, which also can be extracted from the analyzed amplitude data. 

The invention has been described herein with reference to certain examples 

and embodiments. It will, however, be evident that various modifications and 
15 changes may be made to the embodiments described above without departing from 

the scope of the invention as claimed. 
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We claim: 

1 . A method of extracting information about a fluid in a rock comprising: 

a) generating a sequence of magnetic field pulses in the fluid, the 
sequence comprising an initial magnetic field pulse, a first portion that follows the 
initial magnetic field pulse, and a second portion that follows the first portion; 

b) detecting magnetic resonance signals using the second portion of the 
sequence; 

c) modifying the first portion of the sequence, and repeating steps a) and 

b); and 

d) analyzing magnetic resonance signals from a time relative to the initial 
magnetic field pulse and extracting information about the fluid in the rock. 

2. The method of claim 1 , wherein the second portion comprises a series 
of magnetic field pulses separated by a time spacing. 

3. The method of claim 2, wherein the first portion comprises a first 
series of magnetic field pulses separated by a first time spacing. 

4. The method of claim 3, wherein the first time spacing is not less than 
the time spacing of the second portion. 

5. The method of claim 2, wherein the first portion comprises a 
stimulated echo sequence. 
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6. The method of claim 1, wherein the first portion includes at least one 
magnetic field gradient pulse. 

7. The method of claim 1, wherein modifying the first portion comprises 
varying a time spacing between magnetic field pulses. 

8. The method of claim 1, wherein analyzing magnetic resonance signals 
comprises calculating an attenuation factor. 

9. The method of claim 1, wherein analyzing magnetic resonance signals 
comprises determining a relative amplitude of the magnetic resonance signals. 

10. The method of claim 1, wherein extracting information about the fluid 
comprises determining a diffusion coefficient 

1 1 . The method of claim 1 , wherein extracting information about the fluid 
comprises distinguishing between different components of the fluid. 

12. The method of claim 1, wherein extracting information about the fluid 
comprises determining a saturation of the fluid in the rock. 

13. The method of claim 1 , further comprising: 
e) repeating step c) a plurality of times. 
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14. The method of claim 13, wherein extracting information about the 
fluid comprises extracting a two dimensional map of a first parameter indicative of 
the fluid in the rock versus a second parameter indicative of the fluid in the rock. 

15. The method of claim 14, wherein the first parameter is a diffusion 
coefficient. 

16. The method of claim 14, wherein a second parameter is a relaxation 

time. 

17. The method of claim 13, further comprising: 

f) varying the time relative to the initial magnetic field pulse from which 
the magnetic resonance signals are analyzed and repeating steps a) through e). 

18. The method of claim 17, further comprising repeating step f) a 
plurality of times. 



19. The method of claim 17, wherein extracting information about the 
fluid comprises determining a time dependent diffusion coefficient. 

20. The method of claim 17, further comprising determining an indication 
of pore geometry of the rock. 
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21. A method of extracting information about a fluid in a rock comprising: 

a) generating a first sequence of magnetic field pulses in the fluid, the 
first sequence comprising a first portion and a second portion; 

b) detecting magnetic resonance signals using the second portion of the 
first sequence; 

c) generating at least one other sequence of magnetic field pulses in the 
fluid, each other sequence comprising a third portion and a fourth portion, the third 
portion having a different sensitivity to diffusion effects than the first portion of the 
first sequence and the fourth portion having a similar sensitivity to diffusion effects as 
the second portion of the first sequence; 

d) detecting magnetic resonance signals using the fourth portion of each 
other sequence; and 

e) analyzing magnetic resonance signals to separate diffusion effects from 
relaxation effects. 

22. The method of claim 2 1 , wherein the first portion of the first sequence 
is substantially similar to the second portion of the first sequence. 

23. The method of claim 22, wherein the first sequence of magnetic field 
pulses is generated according to: 

90 - [ta^ - 1 80 - t^ - echOiL 
wherein tsbon is a time spacing; echo] is the i* magnetic resonance spin echo; and n is 
the number of spin echoes. 
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24. The method of claim 2 1 , wherein the first sequence comprises a 
modified CPMG sequence. 

25. The method of claim 24, wherein the first sequence is generated 
according to: 

90 - [tjong - 1 80 - - echojjm - [tshoit - 180 - tshoit - echOi«] n - 
wherein tj 0 n g is a first time spacing; echo* is the k 111 magnetic resonance spin echo of 
the first portion; m is the number of spin echoes of the first portion; t^on is a second 
time spacing, the secopd time spacing being shorter than the first time spacing; echo;* 
is the i ,lh magnetic resonance spin echo of the second portion; and n' is the number of 
spin echoes of the second portion. 

26. The method of claim 24, wherein the first sequence is generated 
according to: 

90 - [5 - 90 - (A-8) - 90 - S -echo k ] m - [t ±OIt - 1 80 - t& on - echo r ] n > 
wherein 5 and A are time spacings in the first portion; echoi- is the stimulated 
magnetic resonance spin echo of the first portion; m is the number of stimulated spin 
echoes of the first portion; t^on is a time spacing in the second portion; echor is the i ,t] 
magnetic resonance spin echo of the second portion; and n' is the number of spin 
echoes of the second portion. 

27. The method of claim 21, wherein the other sequence comprises a 
modified CPMG sequence. 



lSDOCID:<WO 0142817A1 IA> 



.WO 01/42817 




S00/32589 



28. The method of claim 27, wherein at least one other sequence is 
generated according to: 

90 - [tiongj - 180 - tiongj ~ echo k j] m j - [tshoit - 180 - tshon - echOi'jL-j 
wherein, for the j* other sequence, t] OJlg j is a third time spacing; echo k j is the k* 
magnetic resonance spin echo of the third portion; (mj) is the number of spin echoes 
of the third portion; tshon is a fourth time spacing, the fourth time spacing being shorter 
than the third time spacing; echorj is the i ,th magnetic resonance spin echo of the 
fourth portion; and (n\j) is the number of spin echoes of the fourth portion. 

29. The method of claim 27, wherein at least one other sequence is 
generated according to: 

90 - [5j - 90 - (Aj-Sj) - 90 - Sj -echo kJ ] mJ - [Wt - 1 80 - t^w - echo r jk j 
wherein, for the other sequence, 5j and Aj are time spacings in the third portion; 
echok j is the k* stimulated magnetic resonance spin echo of the third portion; (mj) is 
the number of stimulated spin echoes of the third portion; tshon is a time spacing in the 
fourth portion; echoi j is the i ,th magnetic resonance spin echo of the fourth portion; 
and (n\j) is the number of spin echoes of the fourth portion. 

30. The method of claim 2 1 , wherein the first portion and the third portion 
include at least one magnetic field gradient pulse. 

3 1 . The method of claim 2 1 , wherein analyzing magnetic resonance 
signals comprises detennining an amplitude of the signal from the other sequence 
relative to an amplitude of the signal from the first sequence. 
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32. The method of claim 21 , wherein analyzing magnetic resonance 
signals comprises analyzing spin echoes from the first sequence that correspond in 
time to spin echoes from the other sequence. 

33. The method of claim 2 1 , wherein analyzing magnetic resonance 
signals comprises analyzing a substantially equal number of spin echoes from the first 
sequence and from the other sequence. 

34. The method of claim 21, further comprising extracting a two- 
dimensional map of two parameters indicative of the fluid in the rock. 

35. The method of claim 21, further comprising generating a static 
magnetic field gradient in the rock. 

36. The method of claim 35, further comprising: 

f) moving the rock and the static magnetic field gradient relative to each 
other; and 

g) repeating steps a) through f) a plurality of times to obtain a profile of 
an attribute through at least part of the rock. 

37. The method of claim 2 1 , wherein analyzing magnetic resonance 
signals comprises correlating spin echoes from the first sequence in time with spin 
echoes from the other sequence and analyzing correlated spin echoes at a given time, 
and further comprising: 
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varying the given time at which correlated spin echoes are analyzed, 



and 



g) 



repeating steps a) through f) a plurality of times. 



38. The method of claim 37, further comprising determining a time 
dependent attribute of the fluid in the rock. 

39. The method of claim 37, further comprising determiuing an indication 
of a pore geometry of the rock. 

40. A method of extracting information about a fluid in rock comprising: 

a) generating a first sequence of magnetic field pulses in the fluid, the 
first sequence comprising a first series of magnetic field pulses with a first time 
spacing; 

b) detecting spin echoes using the first series; 

c) generating a second sequence of magnetic field pulses in the fluid, the 
second sequence comprising a second series of magnetic field pulses with a second 
time spacing and a third series of magnetic field pulses with the first time spacing, the 
second spacing being greater than the first time spacing; 

d) detecting spin echoes using the third series; and 

e) extracting information about the fluid in rock using spin echoes 
detected using the first series and spin echoes detected using the third series. 
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41 . The method of claim 40, wherein extracting information comprises 
analyzing spin echoes detected using the first series that correspond in time to spin 
echoes detected using the third series. 

42. The method of claim 40, wherein extracting information comprises 
analyzing a substantially equal number of spin echoes detected using the first series 
and the third series. 

43. The method of claim 40, wherein extracting information comprises 
determining an amplitude of the spin echoes detected using the third series relative to 
an amplitude of the spin echoes detected using the first series. 

44. The method of claim 40, further comprising: 

f) generating at least one other sequence of magnetic field pulses in the 
fluid, each other sequence comprising a fourth series of magnetic field pulses with a 
time spacing greater than the first time spacing and different from the second time 
spacing and from each other, and a fifth series of magnetic field pulses with the first 
time spacing; and 

g) detecting spin echoes using the fifth series, wherein information about 
the fluid is extracted using the spin echoes detected using the first series, the third 
series, and the fifth series. 
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45. The method of claim 44, wherein extracting information comprises 
analyzing spin echoes detected using the first series that correspond in time to spin 
echoes detected using the fifth series. 

46. The method of claim 44, wherein extracting information comprises 
analyzing a substantially equal number of spin echoes detected using the first series 
and the fifth series. 

47. The method of claim 44, wherein extracting information comprises 
extracting a two-dimensional function of a first parameter and a second parameter 
indicative of the fluid in rock. 

48. The method of claim 47, wherein the first parameter is a diffusion 
coefficient and the second parameter is a relaxation time. 

49. A logging apparatus comprising: 

a logging tool that is moveable through a borehole; and 
a processor that is coupled with the logging tool, the processor being 
programmed with instructions which, when executed by the processor: 



cause the logging tool to: 



i) 



generate a sequence of magnetic field pulses in a region of 



investigation of earth formation surrounding the borehole, the sequence 



comprising an initial magnetic field pulse, a first portion, and a second 



portion; 
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ii) detect magnetic resonance signals produced from the region of 
investigation using the second portion of the sequence; 

iii) modify the first portion of the sequence and repeat steps (i) and 
(ii); and 

cause the processor to: 

iv) analyze magnetic resonance signals from a time relative to the 
initial magnetic field pulse and extract information about the region of 
investigation. 

50. The logging apparatus of claim 49, wherein the second portion 
comprises a series of magnetic field pulses separated by a time spacing. 

51. The logging apparatus of claim 49, wherein the first portion comprises 
a first series of magnetic field pulses separated by a first time spacing. 

52. The logging apparatus of claim 49, wherein the first portion comprises 
a stimulated echo sequence. 

53. The logging apparatus of claim 49, wherein the instructions, when 
executed by the processor, cause the logging tool to:v) repeat step iii) a plurality 
of times. 

54. The logging apparatus of claim 53, wherein the instructions, when 
executed by the processor, cause the processor to extract a two-dimensional function 
of two parameters that characterize the region of investigation. 
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55. The logging apparatus of claim 53, wherein the instructions, when 
executed by the processor, cause the processor to vary the time relative to the initial 
magnetic field pulse from which the magnetic resonance signals are analyzed and 
repeat steps i) through v). 

56. A logging apparatus comprising: 

a logging tool that is moveable through a borehole; and 
a processor that is coupled with the logging tool, the processor being 
programmed with instructions which, when executed by the processor: 
cause the logging tool to: 

i) generate a first sequence of magnetic field pulses in a region of 
investigation of earth formation surrounding the borehole, the first sequence 
comprising a first portion and a second portion; 

ii) detect magnetic resonance spin echoes produced from the 
region of investigation using the second portion of the first sequence; 

iii) generate at least one other sequence of magnetic field pulses in 
the region of investigation, each other sequence comprising a third portion and 
a fourth portion, the third portion having a different sensitivity to diffusion 
effects than the first portion of the first sequence and the fourth portion having 
a similar sensitivity to diffusion effects as the second portion of the first 
sequence; 

iv) detect magnetic resonance spin echoes produced from the 
region of investigation using the fourth portion of each other sequence; and 

cause the processor to: ' ' ' 
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v) analyze detected spin echoes to separate diffusion effects from 
relaxation effects. 



57. The logging apparatus of claim 56, wherein the first portion of the first 
sequence is substantially similar to the second portion of the first sequence. 

58. The logging apparatus of claim 57, wherein the first sequence 
comprises a CPMG sequence. 

59. The logging apparatus of claim 56, wherein the first sequence 
comprises a modified CPMG sequence. 

60. The logging apparatus of claim 56, wherein the other sequence 
comprises a modified CPMG sequence. 

61. The logging apparatus of claim 56, wherein the instructions, when 
executed by the processor, cause the processor to analyze a substantially equal 
number of spin echoes detected using the second portion of the first sequence and 
detected using the fourth portion of the other sequence. 

62. The logging apparatus of claim 56, wherein the instructions, when 
executed by the processor, cause the processor to analyze spin echoes detected using 
the second portion that correspond in time to spin echoes detected using the fourth 
portion. 
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63. The logging apparatus of claim 56, wherein step v) involves analyzing 
spin echoes at a given time, and the instructions, when executed by the processor, 
cause the processor to vary the given time. 

64. A logging apparatus comprising: 

means for generating a sequence of magnetic field pulses in the region of 
investigation, the sequence comprising an initial magnetic field pulse, a first portion 
and a second portion; 

means for detecting magnetic resonance signals using the second portion; 

means for modifying the first portion; and 

means for analyzing magnetic resonance signals from a time relative to the 
initial magnetic field pulse to extract information about the region of investigation. 

65. The logging apparatus of claim 64, further comprising means for 
generating magnetic field gradient pulses in the region of investigation. 

66. A logging apparatus comprising: 

means for generating a first sequence of magnetic field pulses in the region of 
investigation, the first sequence comprising a first portion and a second portion; 

means for detecting magnetic resonance signals from the region of 
investigation using the second portion of the first sequence; 

means for generating at least one other sequence of magnetic field pulses in 
the region of investigation, each other sequence comprising a third portion and a 
fourth portion, the third portion having a different sensitivity to diffusion effects than 
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the first portion of the first sequence and the fourth portion having a similar sensitivity 
to diffusion effects as the second portion of the first sequence; 

means for detecting magnetic resonance signals from the region of 
investigation using the fourth portion of each other sequence; and 

means for analyzing magnetic resonance signals to separate diffusion effects 
from surface and bulk relaxation effects. 

67. The logging apparatus of claim 66, further comprising means for 
generating magnetic field gradient pulses in the region of investigation. 
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